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ABSTRACT 

Radio emission at centimetre wavelengths from highly star-forming galaxies, such as submillimetre galaxies (SMGs), is dominated by 
synchrotron radiation arising from supernova activity. Hence, radio continuum imaging has the potential to determine the spatial extent 
of star formation in these types of galaxies. Using deep, high-resolution (Icr = 2.3 ply beam^'; 0775) centimetre radio-continuum 
observations taken by the Karl G. Jansky Very Large Array (VLA)-COSMOS 3 GHz Large Project, we studied the radio-emitting 
sizes of a flux-limited sample of SMGs in the COSMOS field. The target SMGs were originally discovered in a 1.1 mm continuum 
survey carried out with the AzTEC bolometer, and followed up with higher-resolution interferometric (sub)millimetre continuum 
observations. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them (~ 46 ± 11%) with signal-to-noise 
ratios in the range of S/N = 4.2 - 37.4. Towards four SMGs (AzTEC2, 5, 8, and 11), we detected two separate 3 GHz sources with 
projected separations of ~ 175 - 676, but only in one or two cases (AzTEC2 and 11) they might be physically related. Using two- 
dimensional elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of 0754 ± 0711 for our 18 SMGs 
detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of 4.2 ± 0.9 kpc. No clear 
correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, 
there is a hint of larger radio sizes at z ~ 2.5 - 5 compared to lower redshifts. The sizes we derived are consistent with previous SMG 
sizes measured at 1.4 GHz and in mid-J CO emission, but significantly larger than those seen in the (sub)mm continuum emission 
(typically probing the rest-frame far-infrared with median FWHM sizes of only ~ 1.5 - 2.5 kpc). One possible scenario is that SMGs 
have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-T CO line emission and radio 
continuum emission, and ii) a warmer, compact starburst region giving rise to the high-excitation line emission of CO, which could 
dominate the dust continuum size measurements. Because of the rapid cooling of cosmic-ray electrons in dense starburst galaxies 
(~ 10“* - 10^ yr), the more extended synchrotron radio-emitting size being a result of cosmic-ray diffusion seems unlikely. Instead, 
if SMGs are driven by galaxy mergers - a process where the galactic magnetic fields can be pulled out to larger spatial scales - the 
radio synchrotron emission might arise from more extended magnetised interstellar medium around the starburst region. 

Key words. Galaxies: evolution - Galaxies: formation - Galaxies: starburst - Galaxies: star formation - Radio continuum: galaxies 
- Submillimetre: galaxies 


1. Introduction 

Submillimetre galaxies (SMGs; e.g. ISmail et al. 1997I 
[Hughes et al. 1998] [Barger et al. 1998| ) represent a popula¬ 
tion of distant galaxies where star formation is heavily obscured 
by the dusty interstellar medium (ISM). The star formation 
rates (SFRs) in SMGs lie in the range of ~ 10^-10^ Mq yr“', 
and hence these galaxies stand out as the most intense 
starbursts in the universe (for reviews, see IB lain et al. 20021 
[Casey et al. 2014[l. As the potential precursors to the present-day 


massive elliptical galaxies, SMGs have become one of the pri¬ 
mary targets for understanding galaxy evolution across cosmic 
time (e.g. ISwinbank et al. 20061 IFu et al. 20131 IToft et al. 20141 
[Simpson et al. 201^ . In the context of this evolutionary con¬ 
nection, determining the sizes and size evolution of SMGs is 
crucial. 

Nearly all of the cm-wavelength radio emission from star¬ 
forming galaxies, such as SMGs, is non-thermal synchrotron 
radiation from relativistic electrons accelerated in supernova 
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(SN) remnants produced by the short-lived, high-mass OB-type 
stars (M > 8 Mq; main-sequence lifetime tms ^ 30 Myr). 
Because SNe are tracing the recent/on-going star formation, 
the radio synchrotron emission has the potential to trace the 
spatial scales on which star formation is occurring. This con¬ 
nection between radio emission and star formation is strongly 
supported by the tight infrared (IR)-radio correlation ob¬ 
served in galaxies (e.g. IHelou et al. 19851 IBeck & Golla 198^ 
IXu et al. 19921 ICondonT992l 1^ et al. 20011 IBell 20031 
ITabatabaei et al. 20071 [Murphy et al. 2008[ [Sargent et al. 2010| 
IMoric et al. 20101 IDumas et al. 201 11 1. On the basis of this 
correlation, the IR-emitting region of a star-forming galaxy is 
expected to be comparable in size to that of radio continuum 
emission. However, the most recent studies of the sizes of IR- 
emitting regions of SMGs based on continuum imaging obser¬ 
vations with the Atacama Large Millimetre/submillimetre Array 
(ALMA) show that these are significantly smaller than SMG 
radio sizes presented in the literature ( Simpson et al. 2015a[ 
llkarashi et al. 2015l l. A possible explanation for this discrep¬ 
ancy, as suggested by Simpson et al. (2015a), is cosmic ray 
(CR) diffusion in the galactic magnetic field away from their 
acceleration site, which would render larger radio sizes. To test 
this further here we present a study of radio sizes of SMGs from 
a well selected sample of SMGs in the Cosmic Evolution Survey 
(COSMOS; IScoville et al. 2007t deep field using radio data 
from the Karl G. Jansky Very Large Array (VLA)-COSMOS 
3 GHz Large Project (Icr noise of 2.3 pJy beam \ angular 
resolution 0'.'75; V. Smolcic et al., in prep.). We describe the 
SMG sample and the employed VLA data in detail in Sect. 2. 
The 3 GHz images are presented in Sect. 3, and the analysis 
(size measurements and radio spectral indices) are presented in 
Sect. 4. We compare our results with literature studies in Sect. 5, 
discuss the results in Sect. 6, and summarise the main results of 
the paper in Sect. 7. 

To be consistent with the most recent results from the 
Planck mission dPlanck Collaboration et al. 20 IS) , the cosmol¬ 
ogy adopted in the present work corresponds to the flat ACDM 
universe with the dark energy density Qa = 0.692, total 
(dark-i-luminous baryonic) matter density = 0.308, and a 
Hubble constant of Hq - 67.8 km s * Mpc 


2. Data 

2.1. Source sample 

The target SMGs of the present study - AzTECl-30 - were 
originally discovered in the JCMT/AzTEC 1.1 mm continuum 
survey (18" resolution) towards a COSMOS subheld (0.15 deg^ 
in size) by Scott et al. (2008). The signal-to-noise ratios of these 
SMGs were found to be in the range of S/Ni = 4.0-8.3 (see 
Table 1 in IScott et al. 2008T l. The 15 brightest sources, AzTECl- 
15 (S/Ni imm > 4.6), were imaged (and detected) with the Sub¬ 
millimetre Array (SMA) at 890 ;um (2" resolution) by Younger 
et al. (2007, 2009). More recently, AzTEC16-30 (S/Ni.i„ni = 
4.0 - 4.5) were imaged with the Plateau de Bure Interferom¬ 
eter (PdBI) at 1.3 mm (~ 178 resolution) by Miettinen et al. 
(2015). These interferometric follow-up studies have allowed us 
to accurately determine the position of the actual SMGs giving 
rise to the millimetre continuum emission seen in the single-dish 
AzTEC maps and, in eight cases, to resolve the single-dish emis¬ 
sion into multiple (two to three) components (at ~ 2" resolu¬ 
tion). This way, we can reliably identify the correct 3 GHz coun¬ 
terparts of the target SMGs. We note that even the faintest com¬ 
ponent in our source sample (AzTEC26b) has a 1.3 mm flux den- 
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sity of 0.9 mJy, which corresponds to ~ 4 mJy at the observed- 
frame 850 pm (assuming a dust emissivity index of = 1.5; 
see lMiettinen et al. 2015b . and hence can be considered an SMG 
[cf the classic SMG threshold of S ssoum > 5 mJy refer s to bright 
SMGs (e.g. lHainline et al. 2009IIGonzalez et al. 201 11) 1. We also 
note that none of these SMGs has been detected in X-rays, and 
hence they do not appear to harbour any strong active galactic 
nucleus (AGN) [a typical 3cr upper limit to the flux density in 
the 0.5-2 keV band data of the Chandra COSMOS Legacy Sur¬ 
vey is < 6 X 10“*® erg cm“^ s“* (E. Civano et al., in prep.)]. This 
suggests that the observed radio emission from our SMGs is pre¬ 
dominantly powered by star formation. This is further supported 
by the fact that none of our SMGs were detected with the Very 
Long Baseline Array (VLBA) observations at a high, milliarcsec 
resolution at 1.4 GHz (N. Herrera Ruiz et al., in prep.), yielding 
a 3cr flux density upper limit to S ^ beam * Q 

Our sample of 39 SMGs is listed in Table [T] The coordinates 
given in the table correspond to the (sub)mm peak positions de¬ 
termined in the aforementioned SMA and PdBI studies. Table [1] 
also provides the source redshifts that are based on spectroscopic 
measurements (seven sources), optical to near-infrared (NIR) 
spectral energy distribution fitting (i.e. photometric redshift; 17 
sources), and radio/submm flux density ratios (15 sources). The 
redshift distribution is shown in Eig. [T] We refer to Miettinen et 
al. (2015 and references therein) for further details and discus¬ 
sion on the redshifts of our SMGs. 



Fig. 1. Redshift distribution of the target SMGs divided into three sub¬ 
samples: the red and green histograms show the spectroscopic and pho¬ 
tometric redshifts, respectively, while the blue histogram shows the red¬ 
shift values derived from the radio-to-submm spectral index (see Ta¬ 
ble [TJ. The redshift bins have a width of Az = 0.4. The lower redshift 
limits were placed in the bins corresponding to those values. 


2.2. VLA 3 GHz radio continuum data 

The observations used in the present paper were taken by 
the VLA-COSMOS 3 GHz Large Project (PI: V. Smolcic; 
V. Smolcic et al., in prep.). Details of the observations, data 

' As described in Appendix A, we have detected two 3 GHz sources 
towards AzTECS. The western radio source is associated with our target 
SMG, while the eastern 3 GHz source, physically unrelated to the SMG, 
is also detected at 1.4 GHz with the VLBA. 
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Miettinen et al.: Spatial extent of radio emission in the COSMOS SMGs 


Table 1. Source list. 


Source ID 

Q'2000.0 

^2000.0 

Redshift" 

z reference" 


[h:m:s] 




AzTEC1 

09 59 42.86 

+02 29 38.2 

= 4.3415 

1 

AzTEC2 

10 00 08.05 

+02 26 12.2 

2spec “ 1 • 125 

2 

AzTEC3 

10 00 20.70 

+02 35 20.5 

z,pec = 5.298 

3 

AzTEC4 

09 59 31.72 

+02 30 44.0 

Zpho. = 4.93!“] 

4 

AzTEC5 

10 00 19.75 

+02 32 04.4 

7pho. = 3.05!“] 

4 

AzTECe 

10 00 06.50 

+02 38 37.7 

^radio/submm ^ 3.52 

5 

AzTEC? 

10 00 18.06 

+02 48 30.5 

Zph„, = 2.30 ±0.10 

4 

AzTECS 

09 59 59.34 

+02 34 41.0 

Z.pec = 3.179 

6 

AzTEC9 

09 59 57.25 

+02 27 30.6 

2phot = 

zphof = 2jr]il 

4 

AzTEC10 

09 59 30.76 

+02 40 33.9 

4 

AzTEC ll-N* 

10 00 08.91 

+02 40 09.6 

&pec = 1.599 

7 

AzTECll-S'' 

10 00 08.94 

+02 40 12.3 

&pec = 1.599 

7 

AzTEC12 

10 00 35.29 

+02 43 53.4 

Zpto, = 2.54!“] 

4 

AzTEC 13 

09 59 37.05 

+02 33 20.0 

^radio/submm ^ 4.07 

5 

AzTEC 14-E^ 

10 00 10.03 

+02 30 14.7 

^radio/submm ^ 2.95 

5 

AzTEC 14-W^ 

10 00 09.63 

+02 30 18.0 

Zpho. = 1-30!“? 

_ , _ -3 174-0.29 

2 .phof “ 2 >. 1 / Q '^•7 

4 

AzTEC 15 

10 00 12.89 

+02 34 35.7 

4 

AzTEC16 

09 59 50.069 

+02 44 24.50 

^radio/submm ^ 2.42 

5 

AzTEC 17a 

09 59 39.194 

+02 34 03.83 

Zspec = 0.834 

7 

AzTEC 17b 

09 59 38.904 

+02 34 04.69 

7phot = 

7pho. = 3.00!“« 
Zpho. = 3.20!“] 

5 

AzTEC 18 

09 59 42.607 

+02 35 36.96 

5 

AzTEC 19a 

10 00 28.735 

+02 32 03.84 

5 

AzTEC 19b 

10 00 29.256 

+02 32 09.82 

Zphoi = 1.11 ±0.10 

5 

AzTEC20 

10 00 20.251 

+02 41 21.66 

^radio/submm ^ 2.35 

5 

AzTEC21a 

AzTEC21b 

10 00 02.558 

10 00 02.710 

+02 46 41.74 
+02 46 44.51 

Zpho, = 2.60+“] 
Zppo, = 2.80!“] 

5 

5 

AzTEC21c 

10 00 02.856 

+02 46 40.80 

^radio/submm ^ 1.93 

5 

AzTEC22 

09 59 50.681 

+02 28 19.06 

^radio/submm ^ 3.00 

5 

AzTEC23 

09 59 31.399 

+02 36 04.61 

Zphot = 1.60+“] 

5 

AzTEC24a 

10 00 38.969 

+02 38 33.90 

^radio/submm ^ 2.35 

5 

AzTEC24b 

10 00 39.410 

+02 38 46.97 

^radio/submm ^ 2.28 

5 

AzTEC24c 

AzTEC25^ 

10 00 39.194 

+02 38 54.46 

^radio/submm ^ 3.17 

5 

AzTEC26a 

09 59 59.386 

+02 38 15.36 

Zpho, = 2.50!“] 

5 

AzTEC26b 

09 59 59.657 

+02 38 21.08 

^radio/submm ^ 1.79 

5 

AzTEC27 

1000 39.211 

+02 40 52.18 

^radio/submm ^ 4.17 

5 

AzTEC28 

10 00 04.680 

+02 30 37.30 

^radio/submm >3.11 

5 

AzTEC29a 

10 00 26.351 

+02 37 44.15 

^radio/submm ^ 2.96 

5 

AzTEC29b 

10 00 26.561 

+02 38 05.14 

Zpho. = 1-45+“] 

5 

AzTEC30 

10 00 03.552 

+02 33 00.94 

^radio/submm ^ 2.51 

5 


Notes. The coordinates given in columns (2) and (3) for AzTECl- 
15 refer to the SMA 890 pm peak position ((Younger et al. 2007| 
2009), while those for AzTEC 16-30 are the PdBI 1.3 mm peak 
positions dMiettinen et al. 20151 1. The z^pec, Zphot, and z,-adio/submm 
values are the spectroscopic redshift, optical-NIR photometric red- 
shift, and the redshift derived using the Carilli-Yun redshift indica¬ 
tor dCarilli & Yun 19991 2000). The z references in the last column 
are as follows: 1 = lYun et al. 20151 2 =M. Balokovic et al., in prep.; 
3 = lRiechers et al. 2010l and |Capak et al. 201 1[ 4 = ISmolcic et al. 20T^ 
5 = IMiettinen et al. 20151 6 =D. A. Riechers et al., in prep.; 7 =M. Sal- 
vato et al., in prep. AzTEC 11 was resolved into two 890 pm sources 
(N and S) by Younger et al. (2009). The two components are proba¬ 
bly physically related, i.e. are at the same redshift (see Appendix A). 

AzTEC 14 was resolved into two 890 pm sources (E and W) by 
Younger et al. (2009). The eastern component appears to lie at a higher 
redshift than the western one dSmolcic et al. 2012l l. AzTEC25 was 
not detected in the 1.3 mm PdBI observations dMiettinen et al. 20151 . 


reduction, and imaging can be found in Novak et al. (2015), 
Smolcic et al. (2015a), and V. Smolcic et al. (in prep.). In the 
present paper, we employ - for the first time - the final, full 
3 GHz mosaic imaging of COSMOS (192 pointings in total). 
Briefly, these S-band observations were carried out with the 
VLA of the NRACj^in its A and C configurations (maximum 
baseline of 36.4 km and 3.4 km, respectively) between 2012 
and 2014. The 2 GHz bandwidth (2 basebands of 1 GHz each) 
used was divided into 16 sub-bands/spectral windows (SPWs) 


^ The National Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Eoundation operated under cooperative agreement by As¬ 
sociated Universities, Inc. 


each with a 128 MHz bandwidth. Each SPW was subdivided 
into 64 spectral channels with a width of 2 MHz. The data were 
calibrated using the AIPSLite data reduction pipeline, which is 
an extension of NRAO’s Astronomical Image Processing Sys¬ 
tem (AIPS|^ package (IBourke et al. 20141 K. Mooley et al., in 
prep.), and adapted for the VLA-COSMOS 3 GHz Large Project 
(for details, see V. Smolcic et al., in prep.). Further editing, flag¬ 
ging, and imaging was done using the Common Astronomy Soft¬ 
ware Applications package (CASA[^ IMcMullin et al. 2007b . To 
reduce sidelobes and artefacts in the data, phase solutions ob¬ 
tained from self-calibration with the bright quasar J1024-0052 
were applied on each pointing. Every held was cleaned down to 
5cr, and further cleaned down to 1.5cr using manually defined 
masks around the sources. 

The data used here were imaged using the multi¬ 
scale multi-frequency synthesis (MS-MFS) method 
(IRau & Cornwell 201 111 . Briggs or robust weighting was 
applied to the calibrated visibilities with a robust value of 0.5. 
Considering the aim of the present study (i.e. measuring the 
3 GHz sizes of our SMGs), the main advantage of MS-MFS is 
that the final image resolution is not determined by the lowest 
frequency of the bandwidth used because all the SPWs are 
used in the image deconvolution. A Gaussian u-v tapering was 
applied on each pointing using their own Gaussian beam size 
(Full Width at Half Maximum or FWHM). The final mosaic 
was restored with a circular synthesised beam size (FWHM) of 
0maj = Omin = 0775, where Omni and 0tnin are the major and minor 
axes of the beam. The final Icr root mean square (rms) noise 
level in our maps is typically about 2.3 //Jy beam '. 

To quantify the effect of bandwidth smearing (BWS) in 
our 3 GHz mosaic, we examined the behaviour of the ratio 
of the total integrated source flux density to its peak surface 
brightness as a function of the S/N ratio (e.g. lBondi et al. 20081 
INovak et al. 20151 V. Smolcic et al., in prep.). This comparison 
showed that the effect of BWS in the full 3 GHz mosaic of COS¬ 
MOS is only up to ~ 3%, and no correction for BWS in the 
peak surface brightness is applied in the present study. To further 
examine the importance of BWS, we created images of a sub¬ 
sample of our sources from separate pointings where the source 
distance from the (nearest) phase centre is different. Besides de¬ 
pending on the fractional bandwidth, the magnitude of BWS is 
directly proportional to the angular distance of the source from 
the phase centre. However, no significant radial smearing was 
seen in the aforementioned images, which lends further support 
to negligible BWS. 

3. 3 GHz images and counterpart identification of 
the AzTEC SMGs 

The 3 GHz images towards our SMGs are shown in Fig. We 
note that at the redshifts of our sources, z = 0.834-5.298, we are 
probing rest-frame frequencies of Vrest - 5.5 - 19 GHz (Trest - 
1.6 - 5.5 cm), which are dominated by non-thermal synchrotron 
radiation with the fraction of thermal emission becoming in¬ 
creasingly important at higher frequencies (e.g. ICondon 19921 
(Murphy et al. 20125) 1. The 3 GHz counterparts of our SMGs 

^ http://WWW.aips.nrao.edu/index.shtml 
CASA is developed by an international consortium of scientists 
based at the NRAO, the European Southern Observatory (ESO), the Na¬ 
tional Astronomical Observatory of Japan (NAOJ), the CSIRO Australia 
Telescope National Eacility (CSIRO/ATNF), and the Netherlands Insti¬ 
tute for Radio Astronomy (ASTRON) under the guidance of NRAO. 
See http://casa.nrao.edu 
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were identified by eye inspection of the corresponding images. 
The SMGs AzTECl-9, 11-N and 11-S, 12, 15, 17a, 19a, 21a, 
24b, and 27 are found to be associated with a 3 GHz source (with 
a median offset of 0726; Table |^, i.e. 18/39 or ~ 46% (with a 
Poisson error on counting statistics of ±11%) of our sources are 
3 GHz-emitting SMGs. The S/N ratios of our detected 3 GHz 
sources are in the range of S/N = 4.2 - 37.4, AzTEC7 being the 
most significant detection. We note that the detection S/N ratio 
at 1.1 mm of these 3 GHz-emitting SMGs was found to be in the 
range ofS/Ni = 4.0-8.3 (IScott et al. 2008t . To summarise, 
18 SMGs in our sample are found to be associated with 3 GHz 
emission. A selection of these SMGs, and the additional 3 GHz 
radio sources not analysed further in the present study are dis¬ 
cussed in more detail in Appendices A and B, respectively. The 
SMGs not detected at 3 GHz are discussed in Appendix C. 

4. Analysis 

4. 1. Measuring the size of the radio-emitting region 

We used AIPS package to determine the deconvolved sizes of 
our 3 GHz sources. Two-dimensional elliptical Gaussian fits to 
the image plane data were made using the AIPS task JMFIT. 
The fitting was performed inside a box containing the source, 
and the fit was restricted to the pixel values of > 2.5cr. The 
results are given in Table and illustrated in Fig. To test 
the reliability of our size measurements, we simulated SMGs 
with an assigned size and varying S/N ratio, and fit them in 
the same manner as the real sources. These simulations, de¬ 
scribed in Appendix D, suggest that the sizes provided by JMFIT 
are generally robust within the uncertainties assigned by the fit¬ 
ting task (see Fig. |D.1[ lower panel). As is often done in radio¬ 
continuum surveys, we considered a source to be resolved if its 
deconvolved FWHM size is larger than one-half the synthesised 
beam FWHM (e.g. iMundell et al. 20001 IHo & Ulvestad 20011 
|Urquhart et al. 2009) l. AzTECl, 6 , 15, 19a, 21a, and 24b are re¬ 
solved in both 0niaj and ^min, while AzTEC2 (both components), 
4, 5 (both components), 7, 8 , 11-N, 11-S, and 12 are resolved 
in 0n,aj but unresolved in AzTEC3 and the additional com¬ 
ponent towards AzTEC 8 are unresolved in both axes. The upper 
size limit for unresolved sources was set to one-half the synthe¬ 
sised beam FWHM (< 0738). For AzTEC9, 17a, and 27 only the 
major axis FWHM could be determined by JMFIT, while the fit¬ 
ting task did not provide a value for the minor axis FWHM (the 
output value= 0). In column ( 8 ) in Table]^ we give the projected 
linear FWHM size for those SMGs with known redshift [i.e. not 
just a lower limit to z derived using the Carilli & Yun (1999, 
2000) method]. As mentioned earlier in Sect. 2.2, no correction 
for the negligible BWS in the peak surface brightness or FWHM 
size was applied. 

To calculate the statistical properties of our radio-emission 
size distribution, we applied survival analysis to take the up¬ 
per limits to the size into account. We assumed that the cen¬ 
sored data follow the same distribution as the uncensored val¬ 
ues, and we used the Kaplan-Meier (K-M) method to construct a 
model of the input data [for this purpose, we used the Nondetects 
And Data Analysis for environmental data (NADA: [Helsel 20MT l 
package for R]. However, because more than 50% of our minor 
axis data are censored, the K-M estimator could not be used to 
determine the median value of the minor axis length, and hence 
its value was derived using the Maximum Likelihood Estimator 
(MLE) of the survivor function. The mean, median, standard de¬ 
viation, and 95% confidence interval of the deconvolved FWHM 
sizes are given in Table For example, the median value of 
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the deconvolved ffmaj among the 18 SMGs detected at 3 GHz is 
0754 + 0711 (FWHM), and the median major axis FWHM in lin¬ 
ear units is 4.2 + 0.9 kpc as estimated for our SMGs with known 
redshift (15 sources with either a Zspec or Zphot value available). 
In the subsequent size analysis we will employ the deconvolved 
FWHM of the major axis, because the value of O^aai would set 
the physical extent of a disk-like galaxy, while the minor axis, 
assuming this simplified disk-like geometry, would be given by 
Omia - Oma] X cos(t) [defined so that for a disk viewed face-on 

a = 0°), = 0maj]. 

In Fig.[^ we plot the deconvolved major axis FWHM sizes as 
a function of the 3 GHz flux density (upper panel) and 890 pm 
flux density (lower panel), where for AzTECl-15 was 

taken from Younger et al. (2007, 2009), and for AzTEC17a, 
19a, 21a, 24b, and 27 the value of ^ggo^im was calculated 
from the 1.3 mm flux density by assuming that /J = 1.5 
dMiettinen et al. 20151 1. No statistically significant correlation 
can be seen between the 3 GHz size and the radio or submm flux 
density, but we note that the largest angular major axis FWHM 
sizes are preferentially found among the sources with the lowest 
3 GHz flux densities, although the size uncertainties for those 
sources are the highest. 

4.2. Spectrai index between 1.4 and 3 GHz, and 3 GHz 
brightness temperature 

To further characterise the radio continuum properties of our 
SMGs, we derived their radio spectral index between 1.4 and 
3 GHz (ffij 4 GHz^’ observed-frame 3 GHz brightness tem¬ 

perature (Tb). The 1.4 GHz flux densities were taken from the 
COSMOS VLA Deep Catalogue (ISchinnerer et al. 20101 1 for all 
the sources except AzTECl, 8 , and 11 for which 5 i 4 ghz was 
taken/revised from Younger et al. (2007, 2009); see Col. (2) in 
Table The angular resolution in the 1.4 GHz VLA Deep mo¬ 
saic is 275 (ISchinnerer et al. 2010t . while that of the 1.4 GHz 
VLA-COSMOS Large Project data, used by Younger et al. 
(2007, 2009), is 175 x 174 (ISchinnerer et al. 20071 1. These are 
about 3.3 and 1.9 times poorer than in our 3 GHz mosaic, respec¬ 
tively. This difference was not taken into account, but we used 
the 1.4 GHz peak surface brightness as the corresponding source 
flux density, except for AzTEC 8 and 11, for which Gaussian-fit 
based flux densities from Younger et al. (2009) were used (see 
Table |4]i. The 1.4 and 3 GHz flux densities were then used to 
derive cTj 4 ^^^, where we define the spectral index as oc y". 
The derived spectral indices are listed in Col. (4) in Table the 
quoted errors were propagated from those of the flux densities. 
The 3 GHz Rayleigh-Jeans brightness temperature was calcu¬ 
lated as Tb = c^S vli2ksv^Q.), where c is the speed of light, ^b 
is the Boltzmann constant, and the solid angle subtended by the 
Gaussian source was derived from Q. - n6^^l{4\r\T). The un¬ 
certainties in Tb were derived from those associated with ^sghz 
and the 3 GHz major axis FWHM size [see Col. (3) in Table Q. 
We note that Smolcic et al. (2015a) already derived the values of 
“l 4 gL for AzTECl [-(0.90 ± 0.46)] and AzTEC3 (> -0.09). 
Given the large associated uncertainties, the present spectral in¬ 
dex for AzTECl [-(0.69 + 0.61)] is consistent with the previous 
value, while the lower limit of ^ -0.91 we have derived 

for AzTEC3 is different because of the lower 3 GHz flux density 
determined here. 

In the top panel of Fig. we plot the 3 GHz angular ma¬ 
jor axis FWHM sizes as a function of while the bot¬ 

tom panel shows the 3 GHz Tb values as a function of <+[ 4 ghz- 
Among local luminous and ultraluminous infrared galaxies or 
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Fig. 2. The VLA 3 GHz images towards AzTECl-30 displayed with north up and east left. The greyscale images are shown with power-law 
scaling (except for AzTEC28, where arcsinh scaling is used to better illustrate the intensity scale; a power-law scaling would he completely hlack), 
and the overlaid red contours start from 3cr and increase in steps of Icr except for AzTECS, 7, 8, 11, 12, 22, and 24 where the step is Vz x cr. 
The white dashed contours show the corresponding negative features (starting from -3cr). The plus sign in each panel marks the (suh)mm peak 
position (SMA 890 /rm for AzTECl-15: PdBI 1.3 mm for AzTEC16-30). The black thick ellipse shows the resulting Gaussian fit to the source 
(centred at the peak position, FWHM size, and P.A.). Eor AzTEC2 (both components), 3, 4, 5 (both components), 7, 8 (both components), 11-N, 
11-S, and 12 the size represents an upper limit (see Table|^. Moreover, for AzTEC4 the peak position was not well determined, and for AzTEC9, 
17a, and 27 only the major axis FWHM could be determined by JMFIT. The blue filled circle shows the synthesised beam size (0'.'75 FWHM). 
Note that the areal coverage of the images differs from each other for illustrative purposes; a scale bar indicating the 1" projected length is shown 
in each panel, annotated with the corresponding proper length [kpc] at the quoted SMG redshift (except when only a lower limit to z is available). 


(U)LIRGs, it has been found that smaller sources exhibit a flatter 
radio spectral index jCondon et al. 199T1 Murphy et al. 201 3] l. 
The observed trend of more compact sources exhibiting flatter 
radio spectral indices is an indication of increased free-free ab¬ 


sorption by ionised gas (i.e. free electrons gain energy by ab¬ 
sorbing radio photons during their collisions with ions). No such 
corTelation is obvious in our data, and the lower aj limits 
muddy the interpretation. Also, no obvious trend is found be- 
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ZTEC21 . 
z=2.60 


ZTEC28 
z > 3.11 


AzTEC29b 

z=1.45 


ZTEC30 

z > 2.51 


Fig. 2. continued. 


tween and 40 ^^ as shown in the bottom panel of Fig. 
i.e. sources with a higher do not appear to show spectral 
flattening (cf. Fig. 1 in [Murphy et al. 2013) note their different 
definition of Sy oc v^“). The low values of T^, ranging from 
1.4 + 0.7 K to 75.9 ± 16.6 K, show that the observed 3 GHz radio 
emission from our SMGs is powered by star formation activity 
and no evidence of buried AGN activity is visible in our data 
[AGN have Tg > 3.2 x 10^ K at 8.44 GHz v^est at z ^ 1.8 for 
Vobs = 3 GHzh l^ndon et al. 19911[Murphy et al. 2013| . 


5. Comparison with iiterature 

5.1. Previous size measurements of the COSMOS/AzTEC 
SMGs 

Besides the present work, the sizes of AzTECl, 3, 4, 5, and 8 
have been previously determined at 3 GHz and/or other observed 
frequencies (see Table [^. Below, we discuss the size measure¬ 
ments of these five high-redshift SMGs in more detail. 

AzTECl. We have found that AzTECl is resolved 
(0'.'67/g2o X 0'.'43/q3q) in our 3 GHz image. Smolcic et al. 
(2015a) employed a 3 GHz submosaic of the COSMOS field, 
which was based on 130 hr of observations, and had a Icr rms 
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Table 2. Results of Gaussian fits to the 3 GHz sources. 


Source ID 

fl' 2000 .o“ 

^ 2000.0 ” 

T ^ 

h GHz 

■^SGHz^^ 

S/N 

FWHM size^ 

P.AS 

Offset 


[h:m:s] 


[/iJy beam"* ] 

[/tJy] 


["] 

[kpc] 

[°] 

["] 

AzTECl 

09 59 42.86(±0.003) 

+02 29 38.20(±0.05) 

18.3 ±2.4 

28.3 ±5.2 

8.0 

^•”'- 0.20 ^ 

X 3.0(j j 

118.1(21/ 

0 

AzTECZ'' 

10 00 08.04(±0.003) 

+02 26 12.26(±0.06) 

15.0 ±2.4 

18.9 ±4.7 

6.0 

0.54("/jx < 0.38 

4.6(|;2x < 3.2 

170 9+32.^ 

I /u.>'_32 7 

0.16 


10 00 08.01(±0.007) 

+02 26 10.81(±0.08) 

10.0 ±2.4 

14.0 ±5.0 

4.3 

0.72("/jx < 0.38 



1.51 

AzTEC3 

10 00 20.69(±0.003) 

+02 35 20.37(±0.04) 

19.6 ±2.3 

19.6 ±2.3 

8.5 

<0.38 

< 2.4 


0.20 

AzTEC4 

09 59 31.70(±0.06) 

+02 30 43.96(±0.14) 

11.4 ±2.3 

31.5 ±7.7 

5.3 

1.72(jJ;gx < 0.38 

ll.l+jx < 2.5 

150.6(22 

0.31 

AzTECS'' 

10 00 19.75(±0.001) 

+02 32 04.29(±0.02) 

49.2 ±2.4 

85.8 ±5.8 

21.4 

0.95(“Ujx < 0.38 

7.5(2/x < 3.0 

41.3(j/ 

0.11 


10 00 19.98(±0.003) 

+02 32 09.98(±0.03) 

24.0 ± 2.4 

42.8 ±5.9 

10.7 

1.00+" “x < 0.38 

8.7(| jx < 3.3 

81.3(2/ 

6.56 

AzTEC 6 ^ 

10 00 06.49(±0.005) 

+02 38 37.40(±0.09) 

12.3 ±2.3 

22.4 ±5.8 

5.4 

0.92(“>0.43(|;;“ 

(6.9(j/ X 3.2(“) 

154 3+^^-^ 

0.34 

AzTEC7 

10 00 18.06(±0.001) 

+02 48 30.43(±0.01) 

89.5 ± 2.4 

98.4 ±4.4 

37.4 

0.42(j;;“x < 0.38 

3.5(2jx < 3.2 

29.8(// 

0.07 

AzTECS'' 

09 59 59.33(±0.001) 

+02 34 41.05(±0.02) 

38.8 ±2.4 

49.4 ± 4.8 

16.3 

0.41(Jj"x < 0.38 

3.2('j'*x < 2.9 

46.4(2 

0.16 


09 59 59.51(±0.001) 

+02 34 41.60(±0.01) 

73.5 ±2.4 

73.5 ±2.4 

29.4 

<0.38 

< 3.2 


2.62 

AzTEC9' 

09 59 57.29(±0.002) 

+02 27 30.54(±0.03) 

29.4 ± 2.2 

33.3 ±4.3 

13.0 

ft - - 0 40+0 13 
^maj '-'•^^-0.17 

Q — 'X "J+l l 
^maj ~ '’■“'-1.4 

33.2(2/1 

0.60 

AzTEC 11-N 

10 00 08.90(±0.001) 

+02 40 09.52(±0.01) 

57.0 ±2.3 

67.5 ±4.6 

24.4 

0.40(“^x < 0.38 

3.5(2/x < 3.3 

-5 1 -5+20.5 
•^^••^-20.5 

0.17 

AzTEC 11-S 

10 00 08.94(±0.001) 

+02 40 10.90(±0.01) 

77.5 ± 2.3 

99.6 ±4.8 

33.3 

0.4S(“-“x < 0.38 

4.2(2/x < 3.3 

163 2+*"^-* 

1.40 

AzTEC 12 

10 00 35.30(±0.002) 

+02 43 53.27(±0.02) 

36.5 ± 2.5 

52.5 ±5.2 

14.4 

0.63(";|"x < 0.38 

5.2(2><3.1 

-70 1 +14.8 
'°-^-14.8 

0.20 

AzTEC 15 

10 00 12.95(±0.006) 

+02 34 34.92(±0.10) 

12.2 ±2.4 

27.9 ±6.9 

5.4 

1 j 1 +0.29 ^ A ca+0.26 
^•^^-0.32 ^ ^-^^-0.50 

9.4(/^ X 3.9(22 

146 3+310 

1.19 

AzTEC IVaf 

09 59 39.19(±0.001) 

+02 34 03.58(±0.02) 

34.9 ± 2.3 

40.8 ±4.4 

15.1 

Smaj = 0.45(2“ 

A . _ -I c+ 0.8 

"maj ~ •^•'’-0 9 

34 4+i'^i 

•^^•^-17.1 

0.26 

AzTEC 19a 

10 00 28.72(±0.002) 

+02 32 03.68(±0.03) 

31.6 ±2.2 

45.3 ± 5.0 

14.7 

0.54(2;|i X 0.44(“j 

4.2(“ X 3.4([/ 

174 3+‘*‘*-2 

0.27 

AzTEC21a 

10 00 02.63(±0.01) 

+02 46 42.14(±0.10) 

9.4 ± 2.3 

25.7 ± 7.9 

4.2 

1 -54+0.40 w A ^7+0.32 
-0.44 ^ -0.50 

11.0(2/x5.5(“ 

AC 7+21.8 
^-^•^- 21.8 

1.15 

AzTEC24b" 

10 00 39.28(±0.002) 

+02 38 45.14(±0.03) 

28.0 ± 2.2 

37.6 ±4.7 

12.6 

O 

4^ 

1 + 

P P 

X 

o 

4^ 

1 + 

P P 

(3.8(1/x3.6([/) 

61.3(2 

2.67 

AzTECZ?''-' 

10 00 39.21(±0.004) 

+02 40 52.65(±0.12) 

9.9 ± 2.3 

12.6 ±4.8 

4.3 

n - - 0 Q')+0-34 

Dmaj — *-'•^^-0.41 

(^maj = 

101 - 12 :? 

0.47 


Notes. The meaning of columns is as follows: (1): SMG name; (2) and (3): peak position of the fitted Gaussian; (4): peak surface brightness; 
(5): total flux density provided by the Gaussian fit; (6): S/N ratio as determined from the maximum pixel value with respect to the rms map 
noise; (7) and (8): deconvolved FWHM size x ^min) in arcsec and physical kpc; (9) position angle of the fitted Gaussian measured from 
north through east; (10); projected angular offset from the (sub)mm position.Formal Icr uncertainties in seconds for a 20 oo o and arcseconds 
for ^ 2000.0 returned by JMFIT are given in parentheses. ® The quoted error in /sghz is the Icr rms noise in the map determined inside a 
~ 300 □" box placed near the SMG, and which did not include any 3 GHz sources. The uncertainty in Ssghz represents the formal error 
determined with JMFIT. The uncertainties do not include the absolute calibration uncertainty, The size and RA. uncertainties represent the 
minimum/maximum values as returned by JMFIT. Note that the P.A. is formally defined to range from 0° to 180°, but for example for AzTEC2 
the maximum P.A. value is 203?5, which is equivalent to an angle of 203?5 - 180° = 23?5. The minimum and maximum P.A. values for 
AzTECS and 24b are equal to the nominal value, and hence the quoted uncertainties are equal to zero. Two 3 GHz sources were detected. 
No linear size is reported for the secondary component towards AzTEC2 because of its unknown redshift. For AzTEC6, 24b, and 27 only a 
lower redshift limit is available (see Table[T](, and the linear FWHM size quoted in parentheses was calculated at that lower z limit; these linear 
sizes were not included in the statistical size analysis. ® For AzTEC9, 17a, and 27 only the major axis could be determined by JMFIT (minor 
axis= 0). 


noise level of 4.5 juJy beam”', i.e. about two times higher than 
in the data used here. AzTECl remained unresolved (upper size 
limit was set to < 0'.'7) in the previous map although the angular 
resolution was slightly higher, i.e. 0'.'7 x 076. The 0735 x 0725 
resolution SMA 890 fim observations of AzTECl by Younger et 
al. (2008) showed the source size to be x 0min ~ 073 x 072 
(~ 074 X 073) when modelled as a Gaussian (elliptical disk). 
This Gaussian major axis EWHM at Trest = 167 fim is 2.2 ± 0.6 
times smaller than the size at drest =1.9 cm we have derived (see 
Eig.0. Miettinen et al. (in prep.) used ALMA (PI; A. Karim) 
to observe AzTECl at dobs = 870 /rm (drest = 163 jim) con¬ 
tinuum and an angular resolution of 0730 x 0729. Eitting the 
source in the 870 jim image plane using JMFIT results in the de¬ 
convolved EWHM size of 0739x 0731;|;g Qp which is fairly 
similar to the size derived by Younger et al. (2008) at a compara¬ 
ble wavelength. Based on deep Ultra VISTA observations (~ 078 
resolution at EWHM), Toft et al. (2014) derived an upper limit 
of < 2.6 kpc to the observed-frame NIR size of AzTECl. The 
authors fit two-dimensional Sersic models to the surface bright¬ 
ness distributions, and calculated the effective radius encompass¬ 
ing half the light of the model. This size scale corresponds to a 
Gaussian half width at half maximum (HWHM) size (see Ta¬ 
ble 1 in IToft et al. 2014t . and to be compared with our EWHM 


diameters we multiplied the sizes from Toft et al. (2014) by 2. 
The physical radius reported by Toft et al. (2014) corresponds 
to a diameter of < 0776 in angular units, which suggests that 
the rest-frame UV-optical size of AzTEC 1 could be comparable 
to its EIR size (0maj ~ 073) and/or 1.9 cm radio continuum size 
(0maj = 0767!°;^’). 

AzTEC3. This source is unresolved in our 3 GHz image, sim¬ 
ilarly to that found by Smolcic et al. (2015a) in their 3 GHz im¬ 
age (upper size limit was set to < 077). Riechers et al. (2014) 
used ALMA to observe AzTEC3 at an angular resolution of 
0763 X 0756. In the dobs = 1 mm (drest = 159 ^m) contin¬ 
uum, the deconvolved EWHM size of AzTEC3 was derived to 
be 0740 /qq4 X 0717/gj*, while in the drest = 158 fim [Cii] line 
emission the size was found to be larger, 0763/°°® x 0734/° [°. 
We have derived the drest =16 cm upper EWHM size limit of 
AzTEC3 to be < 0738, which is smaller than the aforementioned 
rest-frame EIR continuum and [Cii] sizes (although the major 
axis EWHM at drest =159 /tm is marginally consistent with our 
upper size limit; see Eig. |^. The observed-frame NIR diameter 
of AzTEC3 derived by Toft et al. (2014) is < 4.8 kpc, i.e. < 0776, 
which is also consistent with our radio emission EWHM size, 
and with the EIR EWHM size from Riechers et al. (2014). 
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Fig. 3. Top: The angular FWHM of the major axis at 3 GHz as a func¬ 
tion of the 3 GHz flux density. Bottom: Same as above but as a function 
of the 890 pm flux density. For AzTEC17a, 19a, 21a, 24b, and 27 the 
value of S 890 ;,m was calculated from S 1,3 mm by assuming that the dust 
emissivity index is j8 = 1.5; these data points are highlighted by red 
filled circles in both panels. The horizontal dashed line marks the me¬ 
dian major axis FWHM of 0754. The upper size limit for AzTEC3 is 
indicated by a downward pointing arrow. 


AzTEC4. For this source, the FWHM size at 3 GHz is de¬ 
termined to be l'.'72^g2gX < 0738. The source appears elon¬ 
gated with a major-to-minor axis ratio of > 3.5, but as shown 
in Fig. the 3 GHz peak position is not well determined by 
JMFIT. The 0786 x 0777 resolution SMA 870 pm observations 
of AzTEC4 by Younger et al. (2010) showed the source size to 
beF„,jXF„i„ = (076+072)x(074+072) [(170+074)x(077+076)] 
when modelled as a Gaussian (elliptical disk). This Gaussian 
major axis FWHM at drest = 147 pm is 2.9^j'2 times smaller 
than the size at drest = 1.7 cm we have derived (see Fig.[^. The 
observed-frame NIR diameter of AzTEC4 derived by Toft et al. 
(2014) is < 5.0 kpc (< 0778), which is consistent with the Gaus¬ 
sian EWHM at rest-frame EIR derived by Younger et al. (2010). 

AzTECS. The 3 GHz EWHM size we have determined for 
this source is 0795j;QQ2X < 0738, i.e. the major axis is re¬ 
solved while the minor axis is unresolved. Miettinen et al. (in 
prep.) used ALMA (PI: A. Karim) to observe AzTEC5 at dobs = 
994 pm (drest - 245 pm) continuum and an angular resolution 
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Fig. 4. Top: The 3 GHz angular major axis FWHM size plotted against 
the radio spectral index between the observed frequencies of 1.4 and 
3 GHz. The arrows pointing right indicate lower limits to while 

the down-pointing arrows show the upper size limits. The data points 
highlighted with red filled circles are for the additional 3 GHz compo¬ 
nents seen towards AzTEC5 (upper data point) and AzTEC8 (lower data 
point); as discussed in Appendix A, those are probably not physically 
related to the SMGs. The horizontal dashed line marks the median major 
axis FWHM size of our SMGs (0754). For reference, the red shaded re¬ 
gion shows the radio spectral index range of -0.8... - 0.7, which is typ¬ 
ical of the non-thermal synchrotron radio emission from star-forming 
galaxies. Within the errors, 11 out of 20 sources (55% ± 17%) shown 
here have a value consistent with this range, and 6 additional 

sources have a lower limit to ctj less than -0.8. Bottom: The 3 GHz 
brightness temperature as a function of a[™Hz- The symbols are as in 
the top panel, except that the arrows pointing up indicate the lower Tg 
limits, and the lower and upper red filled circles are for AzTEC5-N and 
AzTEC8-E, respectively. 


of 0752 X 0730. The source was resolved into two components 
with a projected separation of 0775 (~ 5.86 kpc at the source 
redshift). The northern ALMA component is perfectly coinci¬ 
dent with our 3 GHz source (0703 offset), but we note that the 
3 GHz emission extends towards the southern ALMA EIR com¬ 
ponent, and hence the detected 3 GHz emission encompasses the 
two ALMA-detected components. Eitting the ALMA sources in 
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the image plane using JMFIT results in the deconvolved FWHM 
size of 0''45;^gg2 ^ for the northern component, and 

0"56;^QQg xO"38;^QQg for the southern component. The major axis 
FWHM length of 0''95^Q g 2 3 GHz is comparable to the sum 
of the major axes of the ALMA emission from the two sources 

1"). Toft et al. (2014) used high resolution (FWHM- 0"2) 
data from the Wide Field Camera 3 (WFC3) aboard the Hub¬ 
ble Space Telescope to determine the rest-frame UV/optical size 
of AzTECS. The diameter derived from their reported radius is 
1.0 ± 0.8 kpc (0'.'12 + O'.'IO). The major axis FWHM of AzTEC5 
at 3 GHz is 7.9 ± 6.6 times larger than its UV/optical diameter. 

AzTECS . Eor this source, the EWHM size at 3 GHz is de¬ 
termined to be O'.'TUgj^x < 0738. The major axis is only 
marginally resolved, while the minor axis is unresolved. The 
0'.'86x0'.'55 resolution SMA 870 pm observations of AzTEC8 by 
Younger et al. (2010) showed the source size to be x - 
(076 + 072) X (075 + 073) [(170 + 075) x (074 + 078)] when 
modelled as a Gaussian (elliptical disk). This Gaussian major 
axis EWHM at drest = 208 pm is L5 /q 2 times larger than the 
radio size at drest = 2.4 cm we have derived (see Fig. [^. The 
observed-frame NIR diameter of AzTEC8 derived by Toft et al. 
(2014), < 6.0 kpc (< 0778), is consistent with our 3 GHz size 
and the Gaussian rest-frame FIR size determined by Younger et 
al. (2010). 

For the remaining of our SMGs mostly upper size limits at 
other wavelengths are available. Younger et al. (2007, 2009) 
constrained the observed-frame 890 pm sizes of AzTECl-15 
to < 172 (i.e. the sources were unresolved), with the excep¬ 
tion of AzTECll that was found to be resolved but best mod¬ 
elled as a double point source. Of the 3 GHz-detections among 
AzTEC 16-30, all the other SMGs except AzTEC21a were found 
to be unresolved by Miettinen et al. (2015) in the - 178 res¬ 
olution PdBI 1.3 mm images. A Gaussian fit to AzTEC21a 
yielded a rather poorly constrained deconvolved EWHM of 
(276 + 172) X (073 + 075). Eor a fair comparison with the present 
3 GHz size, we fitted the source using JMFIT, and obtained a 

1.3 mm EWHM size of 2779/“ “ x 0760/° ® (P.A. = 47?8/[‘ °), 
which is comparable to the aforementioned value, and the major 
axis is about 2.1/yg times larger than that at 3 GHz (1734/°-4°). 
AzTEC21a is potentially a blend of smaller (sub)mm-emitting 
sources (IMiettinen et al. 2015b . hence appears more extended at 

1.3 mm than its radio size. We note that the PdBI 1.3 mm emis¬ 
sion of AzTEC27 could not be well modelled by a single Gaus¬ 
sian source model; the major axis EWHM was determined to 
be = 376 (IMiettinen et al. 20151) . Similarly to AzTEC21a, 
AzTEC27 could be a blend of more compact sources. Higher- 
resolution (sub)mm imaging is required to examine the possi¬ 
ble multiplicity of AzTEC21a and 27. The constraints on the 
(sub)mm EWHM sizes of our SMGs are listed in Table [5] a nd 
plotted against their 3 GHz major axis EWHM sizes in Eigd^ 

Besides for AzTECl, 3, 4, 5, and 8, Toft et al. (2014) 
also derived the rest-frame UV/optical sizes for AzTEClO and 
AzTEC15 (see their Table 1). The measurements were based 
on UltraVISTA observations. The diameters were found to be 

1.4 + 0.2 kpc ( 0718 /°° 2 ) ^j. AzTECIO, and 10.0 ± 1.6 kpc 
(1728/°|2) for AzTEC15. The 3 GHz major axis EWHM of 
AzTEC 15 (0maj = ^"^l/g®) is in good agreement with its 
UV/optical extent, while AzTECIO was not detected at 3 GHz. 
We note that heavy obscuration by dust can lead to an apparent 
compact size at the rest-frame UV/optical wavelengths. How¬ 
ever, one would expect the central galactic regions to be more 
extincted compared to the outer portions, which could affect the 
surface brightness profile in such a way that the measured size 


Table 3. The 3 GHz size distribution statistics. 


Parameter 

Value" 

Mean 

0771 ± 0709 (5.5 ± 0.7 kpc) 

Median 

0754 ± 0711 (4.2 ± 0.9 kpc) 

Mean 9^,^ 

0745 ± 0702 (3.3 ± 0.2 kpc) 

Median Smin* 

0735 ± 0704 (2.3 ± 0.4 kpc) 

Standard deviation of 

0739 (2.9 kpc) 

Standard deviation of 9^^^ 

0707 (0.8 kpc) 

95% confidence interval of 

0754 - 0789 (4.0-6.9 kpc) 

95% confidence interval of 

0742 - 0749 (2.9-3.7 kpc) 


Notes. The sample size of the major (minor) axis angular sizes is 18 
(15), while that of the linear sizes is 15 (13), i.e. the number of SMGs 
with either a Zjpec or Zp^ot value available. The median value of 9^^^ 
could not be derived using the K-M estimator. Hence, it was calculated 
using the MLE (assuming lognormal distribution), which is almost iden¬ 
tical to the K-M function. A two-sided 95% confidence interval for 
the mean value computed using the K-M method. 


(e.g. the half-light radius) is larger than in the case of no differ¬ 
ential dust extinction. However, if the extincted outer parts of a 
galaxy fall below the detection limit, the effect might go in the 
opposite direction. 


3.5 


u 3.0 


- 2.5 


2.0 


1.0 


AzTEC17a, 19a, 24b 


AzTECS 


AZTEC5 


0.5 1,0 1.5 2.0 

3 GHz major axis FWHM size [arcsec] 


Fig. 5. The (sub)mm angular major axis FWHM sizes of our SMGs 
[at 890 pm from Younger et al. (2007, 2008, 2009), at 870 pm from 
Younger et al. (2010), or at 1.3 mm from Miettinen et al. (2015)] plot¬ 
ted against their 3 GHz major axis FWHM sizes. For AzTFCl, 3, 4, 
5, 8, and 21a the (sub)mm emission has been (marginally) resolved. 
For AzTFCl, 4, and 8 we plot the Gaussian-fit size from Younger et 
al. (2008, 2010), while the size of AzTFC3 is that determined from 
ALMA observations at 1 mm by Riechers et al. (2014). For AzTFC5, 
the submm size was derived from the 994 pm ALMA image (O. Mi¬ 
ettinen et al., in prep.) using JMFIT (see text for details). Similarly, 
the plotted 1.3 mm size of AzTFC21a was derived by fitting the 
source using JMFIT for a better comparison with the present 3 GHz 
size. The AzTFC27 data point is not shown due to its non-Gaussian 
shape at 1.3 mm (see text for details). Among AzTFC2-15, the upper 
890 pm size limits - marked with arrows pointing down - are < 172 
( [Younger et al. 2007] 2009), while the upper limits to the 1.3 mm sizes 
of AzTFC17a, 19a, and 24b are set to one-half the synthesised beam 
major axis size at FWHM ( IMiettinen et al. 20151 Table 2 therein). The 
red dashed line indicates where the radio and (sub)mm sizes are equal. 
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Table 4. Radio continuum characteristics of the 3 GHz detections. 


Source ID 

S 1.4GHz“ 

[pJy] 

Tb" 

[K] 

^JGHz c 

“ 1.4 GHz 

AzTECl 

48+12 

8.9 ±5.0 

-(0.69 + 0.61) 

AzTEC2 

76+14 

8.8 ±7.4 

-(1.83 + 0.77) 

AzTEC3 

< 30 

> 16.3 

> -0.91 

AzTEC4 

< 36 

1.4 ±0.7 

> -0.91 

AzTEC5 

126+15 

12.9 + 2.1 

-(0.50 ± 0.24) 

AzTEC5-N 

85 + 15 

5.8+ 1.9 

-(0.90 ± 0.45) 

AzTEC6 

< 38 

3.6 ±2.4 

> -1.48 

AzTEC7 

132 ± 22 

75.9+16.6 

-(0.39 + 0.21) 

AzTEC8-W 

102+13 

40.0 ± 26.6 

-(0.95 ± 0.32) 

AzTEC8-E 

160 ± 23 

> 67 

-(1.02 + 0.17) 

AzTEC9 

68+13 

25.0+18.8 

-(1.10 + 0.30) 

AzTECl 1-N 

138 ± 26 

57.4 + 21.9 

-(0.94 ± 0.30) 

AzTECl 1-S 

132 + 26 

58.8+11.4 

-(0.37 ± 0.27) 

AzTECl2 

98+16 

18.0 + 6.0 

-(0.82 ± 0.34) 

AzTECl5 

< 32 

2.6+ 1.5 

> -0.93 

AzTEC17a 

68+13 

23.4+11.6 

-(0.88 ± 0.28) 

AzTECl9a 

78+12 

21.1 ± 10.1 

-(0.71+0.37) 

AzTEC21a 

< 44 

1.9+ 1.4 

> -1.63 

AzTEC24b 

63 + 13 

25.3+18.2 

-(0.68 ± 0.43) 

AzTEC27 

< 39 

2.0+ 1.8 

> -2.63 


Notes. The values of 5 1 . 4 GH 2 were taken from the COSMOS VLA 
Deep Catalogue May 2010 dSchinnerer et al. 20101 except for AzTECl, 
8, and 11. AzTECl exhibits 4cr 1.4 GHz emission, hence is not listed 
in the VLA catalogue, which is comprised of > 5cr sources. The value 
of S 1 , 4 GHz for AzTECl was taken as the peak surface brightness mul¬ 
tiplied by 1.15 to correct for BWS ( [Younger et al. 2007} . The resulting 
flux density of 48 ± 12 ^Jy agrees with the value reported by Younger 
et al. (2007, Table 2 therein). For AzTEC8 and 11 we adopted the 
1.4 GHz flux densities from Younger et al. (2009, Table 2 therein), 
and multiplied them by 1.15 to correct for BWS as noted by the au¬ 
thors. The COSMOS VLA catalogue values of S 1 4 ghz for AzTEC8-W 
and AzTEC8-E are 237 ± 52 pjy and 186 pJy (no error given), respec¬ 
tively. However, AzTEC8-E is a stronger 1.4 GHz source than AzTEC8- 
W ( [Younger et al. 2009| >. For AzTECll the COSMOS VLA catalogue 
gives an integrated flux density value of S 1.4 ghz = 302 ± 45 pJy. Hence, 
we adopted the 5 1 4 ghz values for AzTECl 1-N and 11-S from Younger 
et al. (2009). The 3cr upper limits are reported for the non-detections, 
where Icr == 10 - 14.7 ply bearn^*. ® The value of refers to 
the Rayleigh-Jeans brightness temperature at v^bs = 3 GHz. *7 Radio 
spectral index between the observed-frame frequencies of 1.4 GHz and 
3 GHz. 


5.2. Comparison to SMG sizes from the literature 

In this subsection, we discuss the SMG sizes derived in previ¬ 
ous surveys at different wavelengths. The measured sizes dis¬ 
cussed below are derived using the following four observa¬ 
tional probes: i) radio continuum emission at centimetre wave¬ 
lengths; ii) (sub)mm continuum emission (typically correspond¬ 
ing to rest-frame FIR); Hi) molecular spectral line emission aris¬ 
ing from rotational transitions of CO; and iv) rest-frame optical 
emission tracing the spatial extent of the stellar content. A se¬ 
lection of size distributions derived from the reported data in the 
studies discussed below is shown in Fig. alongside with our 
3 GHz size distribution. 


5.2.1. Ratdio sizes 


A previous work of immediate interest for comparison with 
our results is the MERLIN/VLA 1.4 GHz survey (Icr — 
6 juJy beam^'; ~ 0752 x 0748 resolution) by Biggs & Ivison 
(2008) of the Lockman Hole SMGs (spanning a redshift range 
of Zspec = 1.147 - 2.689). The median deconvolved FWHM size 
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Table 5. Rest-frame FIR/submm and UV/optical sizes of our 3 GHz 
detected SMGs. Besides the angular sizes, the physical sizes are given 
in parentheses when a spectroscopic or photometric redshift is available. 


Source ID 

FIR/submm size^ 

UV/opt. size* 

AzTECl 

0'/3 X 072^ 

< 0776 (< 5.2 kpc) 


(2.1 X l.JkpUy 



074 X 073^ 



(2.8 X 2.1 kpc^y 



A/^YO+O.Ol y A/^-2 1 +0.0l£/ 

-0.01 -0.01 
(2.7/» | X 2.1/g-| kpU/ 


AzTEC2 

< 172 ($ 10.1 kpc) 


AzTEC3 

0'.'40/|;-“ X 0'.'17/“*' 

(2.5+“ X l.l+"-5 kpUy 

< 0776 (< 4.8 kpc) 

AzTEC4 

(0'.'6 ± 0'.'2) X (0'.'4 ± 072/ 

[(3.9 ± 1.3) X (2.6 ± 1.3)kpcT 
(I'.'O ± 074) X (077 ± 076/ 

[(6.5 ± 2.6) X (4.5 ± 3.9) kpU/ 

< 0778 (< 5.0 kpc) 

AzTEC5 

0745/2 “ X0728/S;“« 
(3.5/22x2.2/2-2 kpc?)* 

0712 ±0710 (1.0 ±0.8 kpc) 

AzTECe 

^ 172 


AzTEC7 

S 172 (S 10.1 kpc) 


AzTECS 

(076 ± 072) X (075 ± 073/ 

[(4.6± 1.5)x(3.9±2.3)kpcT 

(170 ± 075) X (074 ± 078/ 

[(7.7 ± 3.9) X (3.1 ± 6.2) kpc^/ 

< 0778 (< 6.0 kpc) 

AzTEC9 

S 172 (S 10.0 kpc) 


AzTECl 1-N 

S 172 (< 10.4 kpc) 


AzTECl 1-S 

S 172 (< 10.4 kpc) 


AzTECl2 

172(5 9.9 kpc) 


AzTECl5 

5 172(5 9.3 kpc) 

1728/222 (10.0 ± 1.6 kpc) 

AzTECl 7a 

S 0793 (S 7.3 kpc) 

AzTECl 9a 

S 0793 (< 7.2 kpc) 


AzTEC21a 

2779/2;“ X 0760/2 22'’ 
(22.9+11 ^ 4.9+22 kpc-)'’ 


AzTEC24b 

S 0788 


AzTEC27 

^maj = 376^ 



Notes. The upper FWHM size limits for AzTEC2-15 at observed- 
frame 890 pm are from Younger et al. (2007, 2009), while those 
for AzTEC17a, 19a, and 24b refer to observed-frame 1.3 mm 
dMiettinen et al. 20151 and represent half of the beam major axis 
FWHM. The diameter at rest-frame UV/optical derived from the ef¬ 
fective radii from Toft et al. (2014) (see text for details). The FWHM 
size derived from SMA 890 pm data by Younger et al. (2008) when 
modelling the source as a Gaussian (upper value) or elliptical disk 
(lower value). The FWHM size measured from the ALMA 870 pm 
image (O. Miettinen et al., in prep.) using JMFIT. A deconvolved 
FWHM size derived through ALMA 1 mm observations by Riechers 
et al. (2014). The FWHM size derived from SMA 870 pm data by 
Younger et al. (2010) when modelling the source as a Gaussian (upper 
value) or elliptical disk (lower value). ® The FWHM size measured 
from the ALMA 994 pm image (O. Miettinen et al., in prep.) using 
JMFIT. The observed-frame 1.3 mm FWHM size of AzTEC21a de¬ 
rived here using JMFIT. The major axis FWHM at 1.3 mm from Mi¬ 
ettinen et al. (2015). 


we derived from their data (their Table 3 of AIPS/JMFIT-derived 
sizes) is (0761 + 0710) x (0731 + 0708)|^This is comparable to 
our median size of (0754 + 0711) x (0735 + 0704) derived from 
2.1 times higher frequency observations. The median linear size 
from Biggs & Ivison (2008), (6.3 + 0.9) x (3.3+ 0.7) kpc^ (scaled 
to the Planck 2015 cosmology), is 1.5 ± 0.4 times larger in the 
major axis than our value of (4.2 ± 0.9) x (2.3 + 0.4) kpc^. If 
we consider only those SMGs in our sample that lie in the red¬ 
shift range studied by Biggs & Ivison (2008), i.e. AzTEC7, 11- 
N, 11-S, and 12, the discrepancy becomes more significant: the 
median linear major axis of these four sources is 3.5 ± 0.5 kpc, 
i.e. 1.8 + 0.4 times smaller than that from Biggs & Ivison (2008). 
However, given the relatively small number of sources in these 

^ The median sizes from other works that we report in this section 
were derived (when relevant) using a survival analysis as described in 
Sect. 4.1. 




















Miettinen et al.: Spatial extent of radio emission in the COSMOS SMGs 


two (sub)satnples [Zspec values are available for eight SMGs in 
the Biggs & Ivison (2008) sample], the latter comparison might 
be susceptible to small number statistics. 

To see how the millimetre flux densities of the Biggs & Ivi¬ 
son (2008) SMGs compare to those of our SMGs, we compiled 
their Bolocam 1.1 mm dLaurent et al. 2005l l. MAMBO 1.2 mm 
(llvison et al. 2005l l. and SCUBA 850 /tm (llvison et al. 2007ll 
flux densities and converted them to 1.1 mm flux densities as¬ 
suming that y 6 = 1.5 when needed. The resulting flux density 
range is S i.i mm = 1 - 4 ];q®- 6 . 0];}'4 mJy, which includes somewhat 
fainter sources than ours with the JCMT/AzTEC 1.1 mm flux 
densities of 5i.imm = 3.3];} ^ - 9 . 3 ]^j -3 mJy dScott et al. 2008b . 
However, the two 1.1 mm flux density ranges are compara¬ 
ble within the uncertainties, and hence the radio size compari¬ 
son is reasonable. Moreover, we found no correlation between 
the Biggs & Ivison (2008) SMGs’ radio sizes and their mil¬ 
limetre flux densities (cf. our Fig. bottom panel). We note 
that Chapman et al. (2004), who studied 12 Hubble Deep Field 
SMGs (z = 1.01 - 2.91) using the MFRLIN/VLA 1.4 GHz 
observations (0'.'2 - 0"3 resolution), found that in most cases 
(67 + 24%) the radio emission is resolved on angular scales of 
~ 1" (~ 8.5 kpc at their median redshift of z = 2.2). The me¬ 
dian diameter (measured above 3cr emission) was reported to be 
0'.'83+0'.'14 (7.0];[ 4 kpc), which is larger than our median 3 GHz 
major axis FWHM by a factor of 1.5 + 0.4, although it was not 
specified whether the median diameter refers to the deconvolved 
FWHM as determined in the present paper [we note that accord¬ 
ing to Biggs & Ivison (2008), the size reported by Chapman et 
al. (2004) is the largest extent within the 3cr contour, hence not 
directly comparable with our FWHM sizes]. The authors con¬ 
cluded that their SMGs are extended starbursts and therefore 
different from local ULlRGs with sub-kpc nuclear starburst re¬ 
gions (e.g. lCondon et al. 1991t . Biggs et al. (2010) used 18 cm 
Very Long Baseline Interferometry (VLBI) radio observations 
at a very high angular resolution of about 30 mas to examine the 
sizes of a sample of six compact SMGs drawn from the Biggs 
& Ivison (2008) sample. Only two of these six SMGs (33 + 6 %) 
were found to host an ultra-compact AGN radio core, and the 
authors concluded that the radio emission from their SMGs is 
mostly arising from star formation rather than from an AGN ac¬ 
tivity. 

5.2.2. (Sub)mm sizes 

Ikarashi et al. (2015) recently derived a size distribution for a 
sample of 13 high-redshift (Zphot ~ 3 - 6 ) SMGs through 1.1 mm 
ALMA observations at ~ 072 resolution. Their SMGs were orig¬ 
inally discovered in the ASTE/AzTEC 1.1 mm observations of 
the SnhwraIXMM-Newton Deep Field (S. Ikarashi et al., in prep.), 
and the reported ALMA 1.1 mm flux densities lie in the range 
of 5i.imm = (1.23 + 0.07) - (3.45 + 0.10) mJy. Compared to 
the JCMT/AzTEC 1.1 mm flux densities of our SMGs, namely 
1.1 mm = 3.3^[ g - 9. 3 ^) 3 mJy (IScott et al. 2008II . the Ikarashi et 
al. (2015) SMGs are fainter, which hampers the direct compari¬ 
son of these two SMG samples. Ikarashi et al. (2015) measured 
the sizes using the uv visibility data directly, and assumed sym¬ 
metric Gaussian profiles. From the values given in their Table 1 
we derived a median FWHM of 0722 + 0704. The linear radii 
reported by the authors are very compact, ~ 0.3 - 1.3 kpc (me¬ 
dian ~ 0.8 kpc), which translate to diameters of ~ 0.6 - 2.6 kpc 
(median ~ 1.6 kpc). These sizes suggest that the high-redshift 
(z > 3) SMGs are associated with a compact starburst region 
(as seen at Trest < 160 - 289 //m), and Ikarashi et al. (2015) 
concluded that the median SFR surface density of their SMGs, 


~ 10^ Mq yr“^ kpc~^, is comparable to that of local merger- 
driven (U)LIRGs and higher than those of low- and high-z (ex¬ 
tended) disk galaxies. 

Simpson et al. (2015a) carried out a high-resolution (0735 x 
0725) 870 yum ALMA survey of a sample of 30 of the brightest 
850 jum-selected SMGs from the SCUBA-2 Cosmology Legacy 
Survey of the UKIDSS Ultra Deep Survey (UDS) held. Their 
target SMGs have 850 yum flux densities of S gso/jm = 8-16 mJy 
(> 4cr; ^i.iinrn = 3.2 - 6.5 mJy if yS = 1.5), and hence are 
mostly comparable to our flux-limited sample (only three of our 
SMGs lie above this flux density range). For a subsample of 23 
SMGs (detected at S/N 870 juni > 10)> they derived a median size 
of 0730 + 0704 (2.4 + 0.2 kp^for the major axis FWHM through 
Gaussian fits in the uv planeNThe authors pointed out that Gaus¬ 
sian fits in the image plane yielded sizes consistent with those 
derived in the uv plane, the median ratio between the two be¬ 
ing FWHM(Mv)/FWHM(image)= 0.9 + 0.2. Given the photo-z 
values of the Simpson et al. (2015a) SMGs, the derived median 
size refers to that at Trest ~ 250 yum. The median angular (linear) 
size at a comparable rest-frame wavelength from Ikarashi et al. 
(2015) is still ~ 36% (~ 50%) smaller than in the Simpson et al. 
(2015a) survey. On the other hand, our median 3 GHz angular 
(linear) major axis FWHM is 1.8 + 0.4 (1.8 + 0.5) times larger 
than the median observed-frame 870 yum size from Simpson et al. 
(2015a). Similarly, Simpson et al. (2015a) concluded that their 
rest-frame FIR sizes are considerably smaller (~ 2 times on av¬ 
erage) than the 1.4 GHz radio-continuum sizes from Biggs & 
Ivison (2008). 

5.2.3. Size of the CO emission 

In Fig. 1^ we also show the sizes of SMGs as derived through 
high-resolution (~ 076) CO spectral line (7 = 3-2 and 7-6) 
observations with the PdBI by Tacconi et al. (2006). We derived 
a median CO-emitting FWHM size of (0740 + 0712) x (0740 + 
0710) or (4.1 + 1.0) X (3.3+ 0.9) kpc^ for the Tacconi et al. (2006) 
SMGs that lie at Zspec = 2.202 - 2.509 and have a reported cir¬ 
cular/elliptical Gaussian-flt (uv plane) FWHM size in their Ta¬ 
ble 1. We note that the target SMGs of Tacconi et al. (2006) have 
SCUBA 850 fim flux densities of S gso/jm = 8.2-10.7 mJy, which 
correspond to 1.1 mm flux densities of S i.imm ~ 3.3 - 4.3 mJy 
(assuming y6 = 1.5). Hence, in terms of S i those SMGs are 
comparable to the faintest sources in our sample (AzTEC21- 
30) where we have only three 3 GHz detections (AzTEC21a, 
24b, and 27). Moreover, three of the Tacconi et al. (2006) 
SMGs appear to be hosting an AGN (SMM J044307+0210, 
J123549+6215, and J123707+6214). Nevertheless, our median 
3 GHz major axis FWHM appears to be comparable to the me¬ 
dian CO-emission major axis FWHM from Tacconi et al. (2006): 
the ratio between the two in angular and linear units is 1.4 + 0.5 
and 1.0+0.3, respectively. The mid- to high-7 CO lines observed 
by Tacconi et al. (2006) are more sensitive to denser and warmer 
molecular gas than lower excitation (7up < 2) lines, and therefore 
the total molecular extent is expected to be larger. Indeed, one of 
the Tacconi et al. (2006) SMGs (J123707+6214) was observed 
in CO(7 =1-0) with the VLA by Riechers et al. (2011), and it 
was found to be more spatially extended compared to that seen in 
CO(3-2) emission. Engel et al. (2010; Table 1 therein) provided 


® Simpson et al. (2015a) do not tabulate the individual source sizes, 
and hence we cannot plot the corresponding dobs = 870 yum size distri¬ 
bution in our Fig.[^ We note that Simpson et al. (2015b) list the angular 
FWHM sizes of these ALMA SMGs in their Table 1, but the source 
redshifts are not tabulated by Simpson et al. (2015a,b). 
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a compilation of different CO rotational transition observations 
towards SMGs, and reported linear HWHM sizes for the SMGs 
as derived using circular Gaussian fits in the uv plane (with two 
exceptions where the quoted size corresponds to the half-light ra¬ 
dius). Their target SMGs are characterised by 850 pm flux den¬ 
sities of SssOfim ^ 5 mJy (5i.itnm ^ 2 mJy if yS = 1.5), and 
this threshold is exceeded by all our SMGs in the JCMT/AzTEC 
1.1 mm survey (IScott et al. 20081 1. From their data we derived a 
median HWHM value of 1.85 + 0.39 kpc, which corresponds to 
a FWHM of 3.70 + 0.78 kpc, consistent with the aforementioned 
size we calculated for the Tacconi et al. (2006) SMGs, and hence 
comparable to our median radio emission size (the median major 
axis FWHM being 4.2 + 0.9 kpc). 

5.2.4. The spatial extent of the stellar emission 

Chen et al. (2015) studied the rest-frame optical sizes of SMGs. 
Based on the HubblefWFC3 observations of a sample of 48 
AFMA-detected z = 1 - 3 SMGs in the Extended Chandra Deep 
Field South (AFESS SMGs), the authors measured a median ef¬ 
fective radius (half-light radius along the semi-major axis within 
which half of the total flux is emitted) of 4.4;|^y^ kpc through 
fitting a Sersic profile to the //igo-band (Tpivot = 1 536.9 nm) sur¬ 
face brightness of each SMG. Simpson et al. (2015a) compared 
their FIR sizes to the optical sizes from Chen et al. (2015), and 
found a large difference of about a factor of four between the 
two (the optical emission being more extended). Interestingly, 
the median radius at the rest-frame UV/optical for the AzTEC 
SMGs from Toft et al. (2014) is only 0.7 kpc (derived using sur¬ 
vival analysis; see our Table for the diameters), but we note 
that most of these sources are very high-redshift SMGs (such as 
AzTEC 1 and AzTEC3), which makes their size determination 
more difficult, and, as mentioned in Sect. 5.1, the measured sizes 
are probably subject to strong dust extinction. 


6. Discussion 

6.1. The spatial extent of SMGs as seen in the radio, dust, 
gas, and stellar emission 

The radio continuum emission, thermal dust emission, and 
molecular spectral line emission can all be linked to the stel¬ 
lar evolution process in a galaxy. Star formation takes place in 
molecular clouds where the gas and dust are well mixed. The 
molecular gas content is best traced by the rotational line emis¬ 
sion of CO. However, different transitions (arising from different 
J levels) have different excitation characteristics, hence are prob¬ 
ing regions of differing physical and chemical properties: the 
high-excitation line emission is arising from denser and warmer 
phase, while low-excitation lines (especially the 7=1-0 transi¬ 
tion) are probing colder, more spatially extended gas reservoirs 
(e.g. llvison et al. 20111IRiechers et al. 2011b . Dust grains absorb 
the UV/optical photons emitted by the young, newly formed stel¬ 
lar population, and then re-emit the absorbed energy in the FIR. 
When the high-mass stars undergo SN explosions, the associ¬ 
ated blast waves and remnant shocks give rise to synchrotron 
radio emission produced by relativistic CRs. This connection is 
believed to lead to the tight FIR-radio correlation (see Sect. 1 and 
references therein). On the basis of this connection, one would 
also expect the FIR- and radio-emission size scales to be similar. 
The galactic-scale outflows driven by the starburst phenomenon 
(SNe, stellar winds, and radiation pressure) are not expected to 
overcome the gravitational potential of the galaxy, hence not dis¬ 
persing the ISM out of the galaxy (this requires a stronger feed- 
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back from the AGN; e.g. lTacconi et al. 2006b . To summarise, the 
radio continuum, rest-frame FIR, and mid- to high-7 CO transi¬ 
tions are all expected to trace regions of active star formation, 
and hence the corresponding spatial extents of their emission are 
expected to be comparable to each other. However, as the size 
comparison in the previous subsection shows, this does not seem 
to be the case for SMGs. 

As discussed in Sect. 5.2, we have found that the 3 GHz 
radio-continuum sizes are comparable to the CO-emission sizes 
from Tacconi et al. (2006) and Engel et al. (2010), but more ex¬ 
tended than the FIR emission seen in other studies, most no¬ 
tably when compared to those from Ikarashi et al. (2015). A 
possible scenario is that SMGs have a two-component ISM: a 
spatially extended gas component, which is traced by low- to 
mid-7 CO line emission and radio continuum emission, and 
a more compact starburst region giving rise to the higher-7 
CO line emission. In the former component, a low dust tem¬ 
perature would lead to a low dust luminosity, while the lat¬ 
ter one - having an elevated dust temperature - could domi¬ 
nate the luminosity-weighted dust continuum size measurements 
(cf. lRiechers et al. 2011b . 

Simpson et al. (2015a) suggested that the larger radio¬ 
continuum size compared to that at rest-frame FIR is the 
result of CR diffusion in the galactic magnetic field (e.g. 
[Murphy et al. 2008) 1. To quantify this, they convolved their me¬ 
dian 870 pm size with an exponential kernel and a scale length 
of 1-2 kpc on the basis of the diffusion length of CR electrons 
in local star-forming galaxies (which is an order of magnitude 
longer than the mean free path of dust-heating UV photons; 
Bicay & Helou 1990) Murphy et al. 2006| 2008). The convolved 


size (FWHM) of 3.8-5.2 kpc they derived is in better agreement 
with the median major axis FWHM of 6.3 + 0.9 kpc from Biggs 
& Ivison (2008; see our Sect. 5.2). However, as pointed out 
by Simpson et al. (2015a), the diffusion scale length of CRs in 
SMGs might be shorter than the aforementioned value because 


of the higher SFR surface density in SMGs (Murphy et al. 2008 
see our Appendix E). 

The rest-frame FIR sizes of AzTEC 1 and AzTEC3 (see 
Sect. 5.1) suggest that they are comparable to their 3 GHz ra¬ 
dio sizes within the uncertainties, but higher-resolution (sub)mm 
continuum imaging of all our SMGs is required to better con¬ 
strain their FIR emission sizes, and to examine whether they rep¬ 
resent the population of very compact SMGs, similarly to those 
from Ikarashi et al. (2015). However, even if the radio size is 
more extended than the FIR emission, the short cooling time of 
CR electrons in starburst galaxies (~ 10^ - 10^ yr) suggests that 
their diffusion through the ISM to spatial scales larger than FIR 
emission is infeasible (see Appendix E for the calculation; the 
diffusion length ranges from only a few tens of pc to ~ 10^ pc). 
Hence, the CR diffusion scenario proposed by Simpson et al. 
(2015a) seems unlikely, and in Sect. 6.2 we will discuss a pos¬ 
sible alternative explanation for a more extended radio emission 
in SMGs. 

A further puzzle is the fact that the rest-frame FIR sizes of 
SMGs appear smaller than the CO-emitting size given the FIR- 
CO correlation found for different types of galaxies at both low- 
and high-z, including SMGs (see e.g. ICarilli & Walter 20131 for 
a review; Fig. 7 therein). The large difference found between 
the rest-frame FIR and optical sizes of SMGs (about a factor 
of four; see our Sect. 5.2.4) led Simpson et al. (2015a) to con¬ 
clude that the spatial extent of ongoing star formation is more 
compact than the spatial distribution of pre-existing stellar pop¬ 
ulation, and that their SMGs might be undergoing a period of 
bulge growth. As pointed out by Chen et al. (2015), if the high- 
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redshift (z > 3) SMGs are progenitors of z ~ 2 compact, quies¬ 
cent galaxies (cQGs; see lToft et al. 2014l l. the high-z SMGs have 
to go through a major transformation to decrease the spatial ex¬ 
tent of the stellar component (and to increase the Sersic index) 
before being quenched. 



Major axis FWHM [kpc] 


Fig. 6. The distrihution of the major axis FWHM sizes of our COSMOS 
SMGs as seen at Vobs = 3 GHz is shown by an open histogram. For 
comparison, the following SMG major axis size distributions are also 
shown: 1.4 GHz sizes from Biggs & Ivison (2008), 1.1 mm sizes from 
Ikarashi et al. (2015), and the extent of CO molecular gas in the SMGs 
studied by Tacconi et al. (2006). The upper size limits were placed in 
the bins corresponding to those values. The vertical dashed lines show 
the corresponding median major axis sizes [4.2 kpc for AzTECl-21a, 
6.3 kpc for the Biggs & Ivison (2008) SMGs, 1.6 kpc for the Ikarashi 
et al. (2015) SMGs, and 4.1 kpc for the CO sizes from Tacconi et al. 
(2006); survival analysis was used to take the upper size limits into ac¬ 
count when calculating the median sizes]. See text for details. 


6.2. Merger-induced extended synchrotron emission 

In the scenario where the FIR size of a galaxy is smaller than 
its radio-emitting region, and where CR electron diffusion - 
due to rapid radiative cooling (~ 10"^ - 10^ yr) - is unlikely 
to be the reason for a more extended radio size (which is po¬ 
tentially the case here; see Appendix E), an alternative inter¬ 
pretation is required. One possibility is that if SMGs are driven 
by mergers (e.g. ITacconi et al. 20081 |Engel et al. 2010) 1, the in¬ 
teracting progenitor disk galaxies can perturb each others mag¬ 
netic fields by pulling them out to larger spatial scales (see 
[Murphy 201 3| l. Hence, a significant amount of non-star forma¬ 
tion related radio emission can arise from the merger system. 
Murphy (2013) concluded that this “taffy”-like merger scenario 
could explain the low EIR/radio ratios and steep high-frequency 
radio spectra of local compact starbursts and those seen in some 
high-z SMGs. In this scenario, mergers are expected to be as¬ 
sociated with stretched magnetic field structures between the 
colliding galaxies, giving rise to synchrotron bridges between 
them and/or tidal tails (ICondon et al. 19931) . The synchrotron- 
emitting relativistic electrons in such bridges might have their 
origin in merger-induced shock acceleration, rather than hav¬ 
ing travelled there from the progenitor galaxies due to the rapid 


cooling time (ILisenfeld & Volk 20101 [Murphy 2013 1 see also 
IDonevski & Prodanovic 2015b . 

The 3 GHz sources investigated here are fairly centrally con¬ 
centrated and no evidence of interaction-induced bridges/tails 
is seen except towards AzTECl, 2, and AzTECll. There is a 
~ 2.6cr 3 GHz feature lying 1'.'5 to the NE of AzTECl, and the 
3 GHz major axis EWHM of AzTECl (0'.'67/° ^^) is larger than 
the sample median major axis EWHM (0'.'54 ± O'.'l 1). AzTEC2 
exhibits an additional 3 GHz source to the SW, which might 
be an indication of a merging pair (or a radio jet). The ad¬ 
ditional source has a major axis EWHM of 0'.'72 /q 44, which 
is also larger than the median of 0754 ± 0711. The two 
3 GHz components seen towards AzTECll share a common 
3cr 3 GHz envelope, but AzTECl 1-N and 11-S both have a 
3 GHz major axis EWHM size smaller than the median value 
(0740/°°’ and 0748/°°4). The 1.4 GHz morphologies of the 
Biggs & Ivison (2008; their Eig. 3) SMGs are generally more 
elongated and clumpy than our sources, which could suggest a 
higher merger fraction among their SMGs, and hence somewhat 
more extended radio emission sizes (see Sect. 5.2.1). However, 
a fair fraction of our target SMGs (~ 36% of the total sam¬ 
ple) show clumpy/disturbed morphologies or evidence of close 
companions at different wavelengths ( Younger et al. 2007) 2009; 
IToft et al. 20141 [Miettinen et al. 2015 1. which could be manifes¬ 

tations of galaxy mergers. 

To conclude, there could be a possible connection between 
merger-driven SMGs and their larger radio-emitting size as com¬ 
pared to EIR emission, as would be expected if the above de¬ 
scribed merger scenario is true. However, the spatial distribution 
of molecular gas, as traced by mid- to high-7 lines, appears to 
be comparable to the Vobs = 3 GHz radio emission size. As de¬ 
scribed in Sect. 6.1, this is to be expected if the observed radio 
size of a galaxy is a direct tracer of its spatial extent of star for¬ 
mation. This would not be consistent with the scenario where the 
CRs emit synchrotron radiation as a result of processes not re¬ 
lated to star formation, such as the aforementioned merger sce¬ 
nario. However, these comparisons between CO and radio me¬ 
dian sizes are, unfortunately, based on measurements obtained 
from different samples and the result can be affected by sub¬ 
tle selection effects. Eor example, the 1.4 GHz radio sizes from 
Biggs & Ivison (2008) are instead larger than the CO sizes from 
Tacconi et al. (2006) (see our Eig.|^, which is qualitatively con¬ 
sistent with the scenario of merger-induced synchrotron emis¬ 
sion. To quantitatively compare the spatial extents of radio emis¬ 
sion and molecular gas component, high-resolution radio and 
CO imaging of the same sample of SMGs is required. 

6.3. Size evoiution as a function of redshift and the effect of 
gaiaxy environment 

In Eig. 1^ we show our deconvolved linear major axis EWHM 
sizes as a function of redshift. No statistically significant corre¬ 
lation can be seen between these two quantities, which is con¬ 
sistent with that found by Simpson et al. (2015a) and Chen et 
al. (2015) at shorter wavelengths. We note, however, that, with 
the exception of AzTEC3 (see below), there is a hint of larger 
radio sizes at z ~ 2.5 - 5 compared to our lower redshift SMGs: 
the z ~ 2.5 - 5 SMGs tend to lie above the median size of 
our sample (4.2 kpc, blue dashed line). Also plotted in Eig. [7] 
are the 1.1 mm EWHM sizes from Ikarashi et al. (2015). These 
authors discussed that the compact sizes of their high-redshift 
(^phot ~ 3 - 6) SMGs support the scenario where they represent 
the precursors of cQGs seen at z ~ 2, which, in turn, are be¬ 
lieved to evolve into the massive ellipticals seen in the present- 
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day (z = 0) universe (IToft et al. 20141) . We note that among the 
Ikarashi et al. (2015) SMG sample, the dobs = 1-1 mm sizes are 
larger at z ~ 3.5-5 than those outside that redshift range, al¬ 
though it should be noted that most of their SMGs in this z range 
have only lower z limits available. As mentioned above, there is 
some resemblance in our data, i.e. the radio sizes appear larger 
at a comparable redshift range of z ~ 2.5 - 5. 

Ikarashi et al. (2015) discussed that if both the radio 
and FIR continuum are tracers of star-forming regions, then 
the z > 3 SMGs are more compact than the lower-redshift 
SMGs typically observed in radio continuum emission (e.g. 
[Biggs & Ivison 2008| l. As shown in Fig. |7j our present VLA 3 
GHz data do not suggest such a trend, and, as mentioned ear¬ 
lier, there is actually a hint of larger radio sizes at z ~ 2.5 - 5 
compared to lower redshifts. However, the highest-redshift SMG 
in our sample, AzTEC3 at z - 5.3, shows the most compact size 
among our sources, consistent with the rest-frame FIR sizes from 
Ikarashi et al. (2015). We note that Capak et al. (2011) found that 
AzTEC3 belongs to a spectroscopically confirmed protocluster 
containing eight galaxies within a 1 arcmin^ area, and therefore 
the environment might also play a role in the galaxy size evo¬ 
lution (see also ISmolcic et al. 2015bl l. However, it is currently 
unclear whether the environmental effects in a galaxy overden¬ 
sity will lead to a more compact or more extented radio-emitting 
size compared to field galaxies. On one side, a protocluster en¬ 
vironment is expected to show an elevated merger rate (e.g. 
IHine et al. 2015b . and, as discussed above, mergers are expected 
to pull the galactic magnetic fields to larger spatial scales, and 
hence lead to a more extended radio synchrotron emission. On 
the other side, the ram and/or thermal pressures of the intra¬ 
cluster medium could compress the ISM of the galaxy, increase 
the magnetic field strength, and hence cause an excess in radio 
emission (consistent with a low IR-radio q parameter of < 2 
for AzTEC3; O. Miettinen et al., in prep.). The aforementioned 
pressure forces can drive shock waves into the ISM, and hence 
accelerate the CR particles ( [Murphy et al. 2009) . Consequently, 
the cooling time and diffusion length-scale of CR electrons can 
decrease (see Appendix E), resulting in a compact radio-emitting 
area. More detailed environmental analysis of SMGs is needed 
to understand this further. 


7. Summary and conclusions 

We have used radio-continuum observations taken by the VLA- 
COSMOS 3 GHz Large Project to study the radio sizes of a sam¬ 
ple of SMGs originally detected with the AzTEC bolometer ar¬ 
ray at 1.1 mm, and all followed up with (sub)mm interferometric 
observations. Our main results are summarised as follows: 

1. Of the total sample of 39 SMGs, 3 GHz emission was de¬ 
tected towards 18 or ~ 46+11% of them (S/N = 4.2 - 37.4). 
Eour sources (AzTEC2, 5, 8, and 11) show two separate 
3 GHz sources. 

2. The median angular radio-emitting size (EWHM) we derived 
is (0754 + O'.T 1) X (0735 + 0704). In linear units, derived for 
the SMGs with known redshift, we obtained a median size of 
(4.2 + 0.9) X (2.3 + 0.4) kpc^. The low brightness temperature 
values ofTs = 1.4+0.7 K to 75.9 + 16.6 K are consistent with 
the radio emission being powered by star formation, rather 
than by an AGN. 

3. We found no obvious cotTelations between the EWHM radio¬ 
size and radio or submm flux density or redshift, which is 
consistent with previous studies at other wavelengths. 
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Fig. 7. The linear major axis FWHM sizes [kpc] of our AzTEC SMGs at 
3 GHz plotted as a function of redshift. The upper size limit of AzTEC3 
is indicated by a downward pointing arrow. Also shown (red points) are 
the 1.1 mm sizes from Ikarashi et al. (2015). The horizontal blue and red 
dashed lines show the corresponding median major axis EWHM values 
of 4.2 kpc and 1.6 kpc, respectively. 


4. We found that our derived 3 GHz sizes are comparable to 
1.4 GHz and CO-emission sizes of SMGs reported in liter¬ 
ature, yet they are ~ 1.7- 2.8 times larger than the median 
rest-frame EIR sizes based on high-resolution ALMA obser¬ 
vations, and reported in literature (see Sect. 5.2 for details). 

5. If both the radio and EIR continuum are tracing the same re¬ 
gions of star formation in a galaxy as expected from the EIR- 
radio correlation, then the differing spatial scales of these 
emissions is puzzling. A possible explanation is that SMGs 
have a two-component ISM: i) an extended gas component 
with a low dust temperature, which gives rise to the low- 
to mid-7 CO line and radio continuum emissions, and ii) 
a warmer, compact starburst region giving rise to the high- 
J rotational line emission of CO, which could dominate 
the dust continuum size measurements. The more extended 
radio-emitting size with respect to the compact EIR-emitting 
size was suggested to be the result of cosmic-ray diffusion by 
Simpson et al. (2015a). However, we have shown here that 
the short electron cooling times of ~ 10“* - 10^ yr in dense 
starburst galaxies do not allow the electrons to spread away 
from their sites of origin to the required spatial scales. Hence, 
it seems more probable that the observed synchrotron emis¬ 
sion partly originates in regions around the active starburst 
region, possibly from extended magnetic fields driven by the 
galaxy merging process. 
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Appendix A: Notes on peculiar 3 GHz sources 


AzTEC2, 5, and 8 each show, in addition to a 3 GHz source 
coinciding with the submm peak, an additional source lying at 
1"51 SW, 6"56 NE, and 2"62 NE from the SMA peak, respec¬ 
tively. The additional 3 GHz feature towards AzTEC2 could, in 
principle, represent a radio-emitting lobe of a jet interacting with 
the surrounding medium or a merger component (projected sep¬ 
aration is 12.7 proper kpc at the redshift of AzTEC2). We note 
that the 1.3 mm emission detected towards AzTEC2 with ALMA 
also shows an additional weak (~ 2.9cr) feature at 2714 to the 
SW of AzTEC2 (Cycle 2 ALMA project 2013.1.00118.S; M. Ar¬ 
avena et al., in prep.), but its peak position lies 0760 away to the 
SW of the 3 GHz feature; this offset is within the large statisti¬ 
cal positional uncertainty of ~ 1" of the weak 1.3 mm feature 
[A0,tat cx (S/N)-']. 

The additional radio source towards AzTEC5 was already 
seen at 1.4 GHz (see Table |^, and it is also visible in the 
Spifzer/IRAC and MIPS images ( Younger et al, 2007| Eig. 1 
therein). This source can be associated with the Herschel- 
detected emission-line galaxy 150.08336-1-02.53619, for which 
a spectroscopic redshift of Zspec(Ha) = 1.42 was reported by 
Roseboom et al. (2012). Also, a 5.4cr detection with ALMA at 
994 pm is obtained towards this source (Cycle 1 ALMA project 
2012.1.00978. S; PI: A. Karim; O. Miettinen et al., in prep.), 
while the ALMA 1.3 mm detection is of ~ 3cr significance 
(M. Aravena et al., in prep.). We note that the most up-to-date 
COSMOS spec-z catalogue gives a lower redshift value of Zspec = 
0.9044, which is based on observations with the Inamori Magel¬ 
lan Areal Camera and Spectrograph (IMACS; M. Salvato et al., 
in prep.); however, the quality flag is 1, i.e. this Zspec is considered 
insecure. Hence, the 3 GHz source NE of AzTEC5 is probably 
a lower-redshift galaxy (AzTEC5 itself has Zphot = 3.05);Q2g; see 
Table T]). Similarly, the two radio sources towards AzTECS were 
already seen at 1.4 GHz ( Younger et al. 2009] Eig. 1 therein; 
see our Table 0: the 1.4 GHz source to the NE of the SMA- 
detected SMG was called AzTEC8.E by Younger et al. (2009), 
and the Spitzer 24 pm emission towards AzTEC8 is coincident 
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with AzTECS.E. This source can be associated with the Chan- 
c/ra/X-ray -detected galaxy CXO C J095959.5-H0234 41 whose 
photo-z is Zphot = 2.420 + 0.060 (ISalvato et al. 201 II I, hence it 
is probably at a lower redshift than the SMG on its western side 
(Zspec = 3.179). This source is also detected with the VLBA at 
1.4 GHz at milliarcsec resolution {S i. 4 ghz = 83.8 piy), showing 
the presence of a radio-emitting AGN (N. Herrera Ruiz et al., in 
prep.; N. Herrera Ruiz, priv. comm.). 

Towards AzTECll, we have detected a double 3 GHz 
source, projectively separated by 1'.'5, and where the south¬ 
ern component is coincident with AzTECll- The northern 
3 GHz source is nearly equidistant from AzTECll-N (1738 
separation) and AzTECll-S (1740 separation). Younger et al. 
(2009) reported that the calibrated visibility data of AzTECll 
show significant structure and are best modelled with a dou¬ 
ble point source (their Table 2). The complex structure of the 
visibility data probably makes the derived source positions to 
be rather uncertain. However, a positional uncertainty of only 
072 in both right ascension and declination for the 890 pm 
peak of AzTECll-N and 11-S was reported by Younger et 
al. (2009; Table 1 therein). The authors also recognised an 
elongated 1.4 GHz source towards AzTECll, where the emis¬ 
sion morphology resembles that seen at 890 pm. Their two- 
component Gaussian fit yielded comparable 1.4 GHz flux den¬ 
sities for the two sources (see our Table 0- The two 3 GHz 
sources seen towards AzTECll share a common radio envelope 
(at the 3cr level), and are probably in the process of merging 
The ALMA 1.3 mm observations at 1736 x 0778 resolution to¬ 
wards AzTECll revealed two SMGs separated by 1745 in pro¬ 
jection (M. Aravena et al., in prep.), and the northern one is well 
coincident with our northern 3 GHz source (only 0711 offset). 

Towards AzTEClS the projected separation between the 
890 pm and 3 GHz positions is fairly large, i.e. 1719. The 
890 pm detection of AzTEClS by Younger et al. (2009) was only 
of moderate significance (4.4cr), and no 1.4 GHz counterpart was 
detected, but it was found to be associated with Spitzer IR emis¬ 
sion. The 890 pm position uncertainty reported by Younger et 
al. (2009) is 073 in right ascension and 072 in declination. The 
1736 X 0778 resolution ALMA 1.3 mm observations towards 
AzTEC15 (M. Aravena et al., in prep.), however, show a perfect 
positional coincidence (0703 offset) with our 3 GHz source of 
5.4cr significance, and hence it is physically related to AzTEC 15. 

The 3 GHz source near AzTEC21a (1715 NE of the PdBI 
position) is our weakest source with a S/N ratio of 4.2. The re¬ 
liability of this 3 GHz source candidate is supported by the fact 
that it is also seen at 1.4 GHz, although the 1.4 GHz source is 
also weak (peak surface brightness of 63 piy beam ' or ~ 3.9cr; 

’ We note that the southern SMA source towards AzTEC 11 was acci¬ 
dentally called AzTEC 11 .N by Younger et al. (2009), while the northern 
component was called AzTEC 1 l.S (see Table 1 in |Younger et al. 2009| l. 
* The spectroscopic redshift of Zspec = 1.599 was measured towards a 
position, which lies 0'.'6 NE of AzTEC 11 -N’s 3 GHz peak position, and 
079 SW from that of AzTEC 11-S (M. Salvato et al., in prep.). The phys¬ 
ical relation of the 3 GHz sources is supported by the very low proba¬ 
bility for a chance association. This can be quantified by calculating the 
Poissonian random probability, P = 1 - , where r is the projected 

angular distance of the two sources, and N is the surface number den¬ 
sity of sources [deg“^] that have flux densities greater than or equal to 
that of the source in question jPownes et al. 198'^IScott & Tout 19891 ). 
Using the VLA-COSMOS 3 GHz Large Project catalogue (V. Smolcic 
et al., in prep.), we estimate that the probability of having a AzTECl 1- 
S type 3 GHz source lying 175 away from that of AzTEC 11-N is only 
P ~ 3.4 X lO^”*. The null hypothesis of chance association is generally 
rejected if P < 5%. 
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see IMiettinen et al. 20151 Pig. A.l therein). Hence, we include 
the 3 GHz source near AzTEC21a in our radio size analysis. 

The 3 GHz source (12.6cr) detected 2767 SW in projec¬ 
tion from AzTEC24b is also detected at 1.4 GHz (Table 0, 
and the 1.4 GHz source was associated with the ASTE/AzTEC 
1.1 mm SMG AzTEC/C48 by Aretxaga et al. (2011), a source 
also detected with Herschel (see IMiettinen et al. 20151) . There is 
also a 1.3 mm-detected ALMA source lying 2763 SW from our 
PdBI detection and 0726 SE from the 3 GHz source in question 
(M. Aravena et al., in prep.). The ALMA detection in particular 
shows that the 3 GHz source is an SMG despite the relatively 
large separation from our PdBI source. We note that there is a 
4.3cr 3 GHz source lying 0747 N of AzTEC27. Miettinen et al. 
(2015) reported the presence of weak 1.4 GHz emission (peak 
intensity of 32.1 piy bearn^* or ~ 2.5cr) towards this position, 
but the nature of this radio emission remains unclear (i.e. noise 
feature or associated with the SMG). The weak 3 GHz source 
associated with AzTEC27 is included in the present radio size 
analysis. The additional 3 GHz radio sources not analysed fur¬ 
ther in the present study are described in Appendix B, and those 
not detected at 3 GHz are discussed in Appendix C. 

Appendix B: Notes on additional 3 GHz sources 

As shown in Pig. 5770 SW from AzTEC20, 8743 NW from 
AzTEC22, and 5751 E from AzTEC23, there is a clearly de¬ 
tected 3 GHz source (10.3cr, 22.8(T, and 5.5cr, respectively). 
Interestingly, these 3 GHz sources are closer to the original 
AzTEC 1.1 mm positions than to the PdBI source candidates 
(07 75 E, 2773 SW, and 4725 NE from A zTEC20, 22, and 23; 
see lScott et al. 200^ IMiettinen et al. 2015t . Morever, the 3 GHz 
sources detected towards the AzTEC20, 22, and 23 fields have 
a Herschel 250 pm detection lying at 2745 SW, 1748 SE, and 
7778 NW, respectively (as based on the cross-correlation with 
the COSMOS SPIRE 250 pm Photometry Catalogue). These ra¬ 
dio sources are not analysed further in the present study. 

Appendix C: Notes on the 3 GHz non-detections’ 
appearances at other wavelengths 

The following 21 SMGs (54 ± 12% of the whole sample) were 
not detected at 3 GHz; AzTEClO, 13, 14-E, 14-W, 16, 17b, 18, 
19b, 20, 21b, 21c, 22, 23, 24a, 24c, 26a, 26b, 28, 29a, 29b, and 
30. 

AzTEClO, 13, 14-E, and 14-W were detected at 890 pm with 
a S/N ratio of 5.3, 4.6, 5.0, and 3.9, respectively, by Younger et 
al. (2009). Only the most significant of these moderate SMA 
890 pm detections, namely AzTEClO, was found to exhibit 
Spitzer IR emission, while none of them was detected at optical 
wavelengths or at 1.4 GHz ( [Younger et al. 2009] l. Hence, their 
non-detection at 3 GHz is not surprising. 

As discussed by Miettinen et al. (2015; Appendix C therein), 
the PdBI 1.3 mm SMG candidates AzTEC16, 17b, 20, 21c, 22, 
24a, 24c, 26b, 29a, and 30 have no multiwavelength counterparts 
and some of them could be spurious. The 1.3 mm S/N ratios of 
these sources were found to be in the range of 4.5-6. Moreover, 
AzTEC29b, although a 7.3cr detection, was found to lie at edge 
of the 1.3 mm map. On the other hand, AzTEC27 and AzTEC28 
(S/Ni 3 njnj = 6 and 5.5, respectively) were among the best 
PdBI detections by Miettinen et al. (2015), but neither of them 
were found to have multiwavelength counterparts; only a trace 
of 1.4 GHz emission (2.5cr) was seen towards AzTEC27 (Ap¬ 
pendix A). AzTEC18, 19b, 21b, 21c, 23, and 26a were detected 
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with S/Ni 3 jj,nj = 4.2 - 9.7, and only AzTEC21c was found to 
not have multiwavelength counterparts. However, none of these 
SMGs was detected at 1.4 GHz. Given the aforementioned prop¬ 
erties, it comes as no surprise that among AzTEC 16-30 there are 
so many 3 GHz non-detections (17 of 22, i.e. 77 + 19%). 


Appendix D: Testing the reiiabiiity of the size 
measurements 

To test the reliability of our EWHM size measurements, we sim¬ 
ulated sources using the CASA (release 4.3.1) Toolkit. We first 
generated mock galaxies with a Gaussian flux distribution, in¬ 
trinsic EWHM size fixed to 0774 x 0745 (i.e. the average de¬ 
convolved EWHM derived for our SMGs where both the ma¬ 
jor and minor axes could be determined or constrained), RA. 
ranging from 0° to 135° in steps of 45°, and flux densities cor¬ 
responding to S/N ratios in the range of 4 to 38 in steps of 
S/N = 2, which cover the observed range of S/N ratios of 
our SMGs (S/N = 4.2 - 37.4). The simulated image was con¬ 
volved to the resolution of 0775 to match the resolution of our 
real data. To obtain a realistic background noise level, the simu¬ 
lated galaxies were added to a 3 GHz map of 1'5 x E5 in size, 
and which was cropped from a source-free region of COSMOS 
(Icr = 2.3 fiiy beam '). The resulting image is shown in the top 
panel in Eig. |D.1| The deconvolved EWHM sizes of the simu¬ 
lated sources were then determined using the AIRS task JMFIT 
as described in Sect. 4.1. 

The bottom panel in Eig. EH shows the ratio of the mea¬ 
sured size to the input size as a function of the S/N ratio. The 
data points are shown separately for the major and minor axes. 
As expected, the size measurement is generally more accurate 
for more significant detections, but within the size uncertain¬ 
ties determined by JMFIT the measured deconvolved sizes are 
in good agreement with the real intrinsic sizes (see the dashed 
line in Eig. |D.1| indicating the one-to-one correspondence). Be¬ 
cause most of our detections are of high significance (median 
S/N = 12.6), these simulations suggest that our size measure¬ 
ments are reliable. 


Appendix E: Calculation of the cosmic-ray electron 
cooling times 

To quantitatively examine the possibility that the observed radio¬ 
continuum sizes of SMGs could be the result of CR diffusion, we 
first calculate the maximum lifetime of the electrons considering 
the radiative energy losses due to synchrotron emission, inverse 
Compton (IC) scattering, bremsstrahlung, and ionisation pro¬ 
cesses. In what follows, we calculate the corresponding cooling 
times using the formulas from Murphy (2009; Sect. 2 therein) to 
which we refer the reader for a more detailed description. 

The redshifts of our 3 GHz detected SMGs range from Zspec = 
0.834 for AzTEC17a to Zspec = 5.298 for AzTEC3. If we assume 
that the critical frequency at which the electron emits most of 
its energy, Vcrit, is equal to Vrest = Vobs(l + z), its value is in the 
range of 5.5 GHz to 18.9 GHz. If the electrons are spiralling 
in a magnetic field whose strength is B ~ 100 /uG (a starburst- 
type B-field; e.g. ILacki & Beck 2013l l. the relationship Vcut = 
1.3 X 10^^(BIjdG)(EIGeY)^ yields CR electron energies of £ ^ 
2 - 3.8 GeV. Eor this case, the synchrotron cooling timescale 
for CR electrons is Tsyn ~ 1.4 x 10®(Vciit/GHz)“'^^(B/yL(G)~^'^^ ~ 
3.2 - 6 X 10^ yr. 

The IC cooling timescale is given by tic ~ 5.7 x 

10^(Vcrit/GHz)-'''2(B/;uG)‘/2(t/rad/10^'2ergcm~3)-i yj.^ 


Simulated SMGs + noise 


O VLA beam 



Fig. D.l. Top: Simulated SMGs added to a noise map to simulate our 
real 3 GHz VLA data. The S/N ratio of the sources increases from left 
to right, top to bottom (S/N -4 - 38), and in each row the P.A. ranges 
from 0° to 135° in steps of 45° (being 0° and 45° for the two bottommost 
objects). The synthesised beam size of 0775 is shown in the bottom left 
comer. Bottom: The ratio of the measured to the simulated input source 
size (deconvolved FWHM) as a function of the S/N ratio. The green 
points show the major axes ratio, while the red points show that between 
the minor axes. The minor axis FWHM for the faintest source (S/N = 4) 
could not be determined by the AIRS task JMFIT. The vertical error bars 
were propagated from the size errors determined by JMFIT. A one-sided 
error bar is shown for those cases where the minimum size could not 
be determined by JMFIT. The horizontal dashed line shows the line of 
equality between the sizes. 


t/rad is the radiation field energy density of the galaxy. To 
estimate the value of t/iad, we adopt a total infrared (8- 
1 000 fim) luminosity range of Lir ~ 10^^ - 10'^ Lq appropri¬ 
ate for SMGs (e.g. Magnelli et al. 201^ ISwinbank et al. 20141 
Ida Cunha et al. 20151 see O. Miettinen et al., in prep., for the 
present SMG sample), and as the characteristic size we use the 
median 3 GHz major axis FWHM size derived here, i.e. 4.2 kpc, 
which corresponds to a circular area of 13.9 kpc^. Using Eq. (4) 
of Murphy et al. (2012a), these values imply t/rad in the range of 
~ 5 X 10“'° - 4.9 X 10“° erg cm“^; for a smaller IR-emitting area 
(Air), the value of t/rad Ajj^* would be higher. The values of VcUt 
and B being as above, we derive rjc ~ 2.6 x 10"' - 4.9 x 10^ yr. In 
the context of IC cooling, it should be noted that our SMG sam¬ 
ple contains high-redshift sources, the most extreme case being 
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AzTEC3 at Zspec = 5.298. At high redshifts, the IC scattering 
between relativistic electrons and the cosmic microwave back¬ 
ground (CMB) - boosting the CMB photon energy - becomes 
more important compared to the low-z universe. The reason for 
this is that the energy density of the CMB increases steeply with 
redshift, namely t/cMB (1 + For instance, at Zspec = 5.298, 
the CMB photon bath has about 140 times higher energy density 
than that at the lowest-redshift SMG in our sample (AzTEC17a 
at Zspec = 0.834). This means that besides the intense IR radia¬ 
tion field in a starburst, 1C losses off the CMB photons have the 
potential to increase the cooling rate of the CR electrons (e.g. 
Lacki & Thompson 20101. 

The bremsstrahlung lifetime is Tbrem ~ 8.6 x 

10^(nH/citi ^ yr, where uh is the hydrogen number density 
of the ISM. Assuming a typical average density range of 
uh = 10^ - 10^ cm“^, we obtain Tbrem ~ 8.6 x 10"^ - 8.6 x 10^ yrj^ 

The timescale for ionisation losses can be written as 

Tion ~ 3.6 X 10'°(Vcrit/GHz)'/2(B/^G)-'/2(nH/cm-3)-i X 

[3/2 X ln(vciit/B) -1-49] * yr, which under our assumptions lies 
in the range of Tion ~ 1.9 x 10^ - 3.4 x 10^ yr. 

Finally, due to the combined energy losses from the afore¬ 
mentioned processes, the effective cooling lifetime for CR elec¬ 
trons is given by 


the wind (e.g. lVolk et al. 199^ . The effect of losses due to elec¬ 
tron advection would further shorten the diffusion length-scale. 


fcool “ 1 1 • (E.l) 

f svn + 'Zir + "Zk -I- T. * 

IL brem ion 

The individual timescales calculated above yield TcooI ~ 1.7 x 
10"^ - 1.9 X 10^ yr. In the case of random-walk diffusion, the 
electrons’ escape scale-length is given by /esc = (f^ETesc)'^^, and 
when the diffusion coefficient De is in the energy-dependent 
regime (£ > 1 GeV), the escape length is given by /esc ~ 
7.1 X 10-4(Tesc/yr)'/2(Vcrit/GHz)i/2(B/juG)-'/2 kpc. During the 
above derived cooling time (tcsc = Tcooi) the electrons can 
travel only /esc ~ 22 - 135 pc. Hence, we conclude that if the 
FIR/star-forming sizes of our SMGs are as compact as those 
from Simpson et al. (2015a; 2.4 + 0.2 kpc in median FWHM) 
and Ikarashi et al. (2015; ~ 1.6 kpc in median FWHM), it seems 
unlikely that CR electrons would have had time to propagate 
from their sites of origin to the large distances where we ob¬ 
serve the 3 GHz emission [the median major axis FWHM size 
being 4.2 ± 0.9 kpc]. In the above analysis we did not add 
the effect of the IC scattering off the CMB radiation, although 
at high redshift it can shorten the electron lifetime and diffu¬ 
sion length scale even more. However, apart from AzTEC3, our 
Fig. 0 does not show evidence of smaller radio sizes at higher 
redshifts as expected if the electrons have less time to travel 
away from their site of origin. This is possibly a manifestation 
of the fact that in starburst galaxies, at whatever redshift they 
might be, the local stellar radiation field is intense (cf. the above 
estimate), and hence the CR electrons can suffer from strong 
IC losses from stellar light besides/instead from the CMB (e.g. 
ILisenfeld et al. 199^ [Lacki & Thompson 2010| l. Moreover, we 
have ignored the fact that if the galaxy is associated with a 
galactic-scale wind, the CR particles in the wind adiabatically 
lose momentum and energy on the course of the expansion of 


^ We note that under the assumption of a magnetic flux freezing, 
i.e. B oc (e.g. ICrutcher 1999b . our adopted field strength of 
100 pG would imply a density of «h - 317 cm“^ if S = 10 pG at 
«H = 1 cm“^ as is typically observed in normal star-forming galaxies 
(cf. [Murphy et al. 2009^ . Our adopted range of densities brackets this 
value of wh- 
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